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ABSTRACT 

We present the results of a Nobeyama 45-ni H2O maser and NH3 survey of all 94 northern GLIMPSE 
Extended Green Objects (EGOs), a sample of massive young stellar objects (MYSOs) identified based 
on their extended 4.5 /im emission. We observed the NH3(1,1), (2,2), and (3,3) inversion lines, and 
detect emission towards 97% , 63%, and 46% of our sample, respectively (median rms ~ 50 mK). 
The II2O maser detection rate is 68% (median rms ^ 0.11 Jy). The derived H2O maser and clump- 
scale gas properties are consistent with the identification of EGOs as young MYSOs. To explore the 
degree of variation among EGOs, we analyze subsamples defined based on MIR properties or maser 
associations. H2O masers and warm dense gas, as indicated by emission in the higher-excitation 
NII3 transitions, are most frequently detected towards EGOs also associated with both Class I and 
II CH3OH masers. 95% (81%) of such EGOs are detected in H2O (NH3(3,3)), compared to only 
33% (7%) of EGOs without either CH3OH maser type. As populations, EGOs associated with Class 
I and/or II CH3OH masers have significantly higher NH3 linewidths, column densities, and kinetic 
temperatures than EGOs undetected in CH3OH maser surveys. However, we find no evidence for 
statistically significant differences in H2O maser properties (such as maser luminosity) among any 
EGO subsamples. Combining our data with the 1.1 mm continuum Bolocam Galactic Plane Survey, 
we find no correlation between isotropic H2O maser luminosity and clump number density. H2O maser 
luminosity is weakly correlated with clump (gas) temperature and clump mass. 

Keywords: infrared: ISM — ISM:jets and outflows — ISM: molecules — masers — stars: formation 



1. INTRODUCTION 

The early stages of massive star formation remain 
poorly understood, due in part to the difficulty of iden- 
tifying young massive young stellar objects (MYSOs) 
that are actively accreting and driving outflows. Large- 
scale Spitzer Space Telescope surveys of the Galac- 
tic Plane have recently yielded a promising new sam- 
ple of candidates : Ex t ended Green Objects (EGOs; 
ICveanowski et al.l 120081 120091 ) . selected based on ex- 
tended 4.5 fxm emission, and named for the common 
coding of three-c olor InfraRed Array Camera (IRAC; 
iFazio et al.l [200l images (RGB: 8.0, 4.5, 3.6 /zm). 
Modeling, mid-infrared (MIR) spectroscopy, and nar- 
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rowband near-infrared (NIR) imaging have shown that 
shock-excited molecular line emission, predominantly 
from II2, can dominate the 4.5 / im broadband flux in 
active proto s tellar outflows (e.g . ISmith fc RosenI 120051: 
Smith et all [20061: jPavis et all [20071: lYbarra fc Ladal 



20091 : lYbarra et al.l 120101: IDe Buizer fc Vaccal I20m 
While all the IRAC filters include H2 lines, only the 4.5 
fim band lacks Polycyclic Aromatic Hydrocarbon (PAH) 
emission features (e.g. Fig. 1 of [Reach et al. I [2006) . 
which are readily excited in massive star forming regions 
(MSFRs). Morphologically distinct extended 4.5 ^m 
emission is thus a common feature of well-known MSFRs 
(e.g. DR21, S255N, NGC6334I(N), G34.4-F0.23, IRAS 
18566-H0408: iDavis et al . 2007; Cvganow ski"eral] \200% 
Hunter et al.l [20061: [Shepherd et al. 20071: [Arava et al.l 
2007D, and a means of identifying candidate MYSOs with 
active outflows. 

Cyganowski et al. (2008, hereafter IC08I ) cataloged 
over 300 EGOs in the Galactic Legacy Infrared Mid- 
Plan e Survey Extraordin aire survey area (GLIMPSE- 
I; C hurchwell et al.l [2009D . At the time, the only data 
available for most EGO s were IR surveys. Using the 
GLIMPSE images, divided cataloged EGOs into 
"likely" and "possible" outflow candidates based on the 
morphology and angular extent of t heir e xtended excess 
4.5 /im emission. As detailed by [C08[ . two phenom- 
ena in the IRAC images have the potential to be con- 
fused with moderately extended 4.5 /im emission: multi- 
ple nearby point sources and image artifacts near bright 
IRAC sources. To categorize the [col EGOs, two ob- 
servers independently reviewed three-color IRAC images: 
if either observer thought the MIR morphology could be 
attributable to one of these phenomena, the EGO was 
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considered a "pos sible" outflow candidate. Of the 302 
EGOs in the [COl catalog, 133 (44%) were classified as 
"likely" outflow candidates, 165 (55%) as "possible" out- 
flow candidates, and 4 (1%) as "outflow-only" sources 
(in which the extended outflow emission could be readily 
separated from the central source). ,C08. also tabulated 
whether each EGO was or was not associated with an 
Infrared Dark Cloud (IRDC) visible against the diffuse 8 
/zm background. A majority (67%) of GLIMPSE EGOs 
are associated with IRDCs, which arc thought to be sites 
of the earliest stages of massive star and cluster forma- 
tion Ce.g. iRathb ornc ct al. 2006; Rathbo rne et al.ll20d'^ 
[Chambers et al., ,2009.; Wang et al.. .2011.) . A somewhat 
higher fraction of EGO "likely" outflow candidates are 
found in IRDCs: 71% compared to 64% of "possible" 
outflow candidates (|C08l ). The GLIMPSE survey is too 
shallow to detect distant low-mass outflows; based pri- 
marily on the MIR data, ,C08 argued that GLIMPSE 
EGOs were likely outflow-driving massive YSOs. 

Testing this hypothesis required correlating extended 
4.5 /xm emission with other massive star formation trac- 
ers at high angular resolution. Interferometric studies 
at cm-mm wavelengths have provided much of the key 
evidence to date that EGOs are indeed young, massive 
YSOs driving active outflows. The first strong evidence 
was remarkably high detection rates for two diagnostic 
types of CH3OH masers in sensitive, high angular resolu- 
tion Very Large Array (VLA) surveys (Cyganowski et al. 
2009, hereafter C09): 6.7 GHz Class II and 44 GHz Class 
I CH3OH masers. Radiatively pumped by IR emission 
from warm dust. Cl ass II CH3OH masers are excited near 
the (p roto)star (e.g. lCragg et al.|[2005HCvganowski et al.l 
120091 and references therein), and recent work suggests 
that the luminosities and relative strengths of different 
Class II transitions change as the central source evolves 
(e.g. lEllingsen etHI 120111: iBreen et al.l 120 111 and refer- 
ences therein). The 6.7 GHz transition is the strongest 
and most common Class II CH3OH maser; importantly, 
numerous searches have shown that these masers are 
not f ound towards low- i nass YSOs (e. g. Minicr ct al. 
200I iBourke et al.l [20051 iXu et al.l[200l fPandian et al.l 
20081 ). Collisionally excited in the presence of weak 



shocks, Class I CH3OH masers are generally associ- 
ated with molecula r outflows and outflow/cloud interac- 
tions (e.g. Plambcc FfcMentenl [19901: iKurtz et al.l 120041: 
iVoronkovet al . 2006) , though recent work suggests Class 
I masers may also b e excited by shocks d riven by ex- 
panding HII regions (jVoronkov et al.l[2010[ ). As a result 
of their association with outflows. Class I CH3OH masers 
are more spatially distributed than Class II masers, and 
may be found m any t ens of arcseconds from the driving 
(proto)star fe.g. lC09( ). 

^C0§ detected 6.7 GHz CH3OH masers towards >64% 
of their 28 EGO targets, and 44 GHz CH3OH masers to- 
wards ~90% of the subset searched for Class I emission 
(19 EGOs, 18 with 6.7 GHz CH3OH masers). Their full 
sample of 28 EGOs was chosen to be visible from the 
northern hemisphere and to span a range in MIR prop- 
erties including presence/absence of 8 and 24 fim coun- 
terparts, morphology, IRDC association and angular ex- 
tent of 4.5 /im emission. The 19 sources observed with 
the VLA at 44 GHz were all "likely" outflow candidates 
and, in essence, a 6.7 GHz CH 3OH maser-selected sub- 
sample (for further details see IC09( ) . Subsequent high- 



resolution mm-A observations of two of the IC09I EGOs 
revealed high- velocity bipolar molecular outflows coinci- 
dent with the 4.5 /xm lobes, driven by compact millime- 
ter continuum cores that exhibit hot core line emission 
(iCvganow ski et 3111201 laL hereafter CI la). Recently, ex- 
ceptionally deep VLA 3. 6 an d 1.3 cm continuum obser- 
vations of a sample of 14 IC09 EGOs have shown that the 
vast majority of the tar gets (12/14) are not ultr acom- 
pact (UC) HII regions (jCvganowski et alll201lU here- 
after Cllb). Most (8/14) are undetected at both 3.6 and 
1.3 cm (cr ^ 30 and 250 fiJy beam"^, respectively); four 
sources are associated with weak (<1 mJy) cm-A emis- 
sion consistent with hypercompact (HC) HII regions or 
ionized winds or jets. Based on their cm survey results 
and complementary multiwavelength data, iCllbl argued 
that these EGOs represent an early stage of massive star 
formation, before photoionizing feedback from the cen- 
tral MYSO becomes significant. 

Detailed, high-resolution followup studies have, of ne- 
cessity, been limited to relative ly sm all EGO subsamples, 
and have generall y foc used on lC08l "likely" outfiow can- 
didates (see also >C0§) . Assessing the variation within 
the C08i catalog and the significance of their MIR clas- 
sifications requires large, uniform surveys in tracers of 
dense gas and star formation activi ty. Few su c h sur - 
veys have been conducted to date. iChen et all ()2010l ) 
searched 88 (of 94) northern {S > -20°) EGOs for 
3 mm HCO+, ^^CO, "CO, and C^^^O emission, with 
the primary goal of detecting infall signatures. They 
found a larger "blue excess" towards EGOs associated 
with IRDCs compared to those not associated with 
IRDCs, and towards "possible" compared to "likely" out- 
fiow candidates; however, the interpretation of these re- 
sults was complicated by the likelihood that multiple 
sources/dynamical phenomena we re present within their 
large (~60-80") beam. Recently, iHe et all (|2012| ) con- 
ducted a 1 mm line survey, covering ^251.5-252.5 GHz 
and ^260.2-261 2 GHz, t o wards 89 northern EGOs (res- 
olution -29"). IHe et al.l (j2012D focus on hncwidth and 
line luminosity correlations, however, and do not an- 
alyze EGO subsamples. Chen et al. (2011, hereafter 
CEll) searched for 95 GHz Class I CH3OH masers to- 
wards 192 northern and southern EGOs (of 302 total) 
with the MOPRA telescope (6Ifwhp--36", 3cr~1.6 Jy). 
They found a higher 95 GHz CH3OH mase r det ection 
rate towards "likely" than towards "possible" IC08I EGOs 
(62% and 49%, respectively), and very similar detec- 
tion rates towards EGOs associated/not associated with 
IRDCs (55%/53%). Their Class I CH3OH maser de- 
tection rate is also much higher towards EGOs associ- 
ated with Class II CH3OH masers (80%) than towards 
those without (38%), consistent with the very high Class 
I maser detection rate of IC09I . 

Like Class I C H3OH masers, H2O masers are collision- 
ally pumped (e.g. lElitzur et al.ll989D and associated with 
protostellar outfiows; notoriously variable, H2O masers 
also often exhibit high-velocity emission features, offset 
by 30 km s~^ or more from the systemic velocity (e.g. 
IBreen et al.llMObI : ICaswU fc Breenl[20T?l . While Class 
I CH3OH masers are excited under moderate condi tions 
(T-80 K, n(H2)~105-106 cm'^, e.g. ILeurinill200l and 
associated with outffow-cloud interfaces, H2O masers re- 
quire more extreme conditions (T— 400 K, n(H2)— 10*- 
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10^° cm-'^. lElitzur et al.|[T989l) and are thought to orig- 
inate behind fast shocks in the inner regions of the out- 
flow base. Numerous correlations have been reported 
between the properties of H2O masers and those of the 
driving source or surrounding clump, including recent ev- 
iden ce that Lff^p oc Lhoi over many orders o f magnitude 
(e.g. lUrquhart et"aLl I20T1I: iBae et al.l [20Tl[) . This sug- 
gests that H2O masers may be used to investigate the 
properties of their driving sources, at least in a statisti- 
cal sense for different subsamples-a possibility of inter- 
est for EGOs, since their bolometric luminosities are in 
most cases poorly constrained by available data (see also 

IcTIH) . 

Large H2O maser and NH3 surveys with single- 
dish telescopes have long been recognized as power- 
ful tools f or characterizing mass i ve star forming re- 
gions (e.g. ChurchwcU elaij 11990|; lAnglada et all 119961: 
ISridharan et al. _2 002,), and continue to be applied to 
new sa mples (e.g. lUrquhart et al.l 1201 It iDunham et al.l 
2011bf ). NH3 traces high-density gas (~10^ cm~^, e.g. 



Evansl[T999 : Stabler & Palla 2005) , and provides a wealth 



of information about clump kinematics and physical 
properties; notably, it is an excellent "thermometer." 
This paper presents the results of a H2O maser and NH3 
surv ey of the 94 northern {5 > —20°) EGOs from the 
IC08I catalog with the Nobeyama Radio Observatory 45-m 
telescope. The motivation for this survey was to charac- 
terize the properties of the COS EGO sample as a whole, 
the main goals being to evaluate the significance of the 
MIR classifications from ] COS and to place EGOs in the 
context of other large MYSO samples. We also compare 
the H2O maser and NH3 properties of EGO subsamples 
associated with Class I and/or II CH3OH masers and ex- 
plore correlations between H2O maser and clump prop- 
erties. Evolutionary interpretations have been suggested 
for both CH3OH masers and HpO maser properties (e.g . 
Ellingsen|[2006t lEllingsen et al.|[2007l : IBreen et al.ll2010al: 



Breen fc Ellingsenll2011i r and our survey, in conjunction 



with the 1.1 mm Bol ocam Galactic Plane Survey (BGPS; 
lAguirre et al.l 120111 : [Rosolowsky et al. 2010). provides 
the necessary data to test these scenarios. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Nobeyama 45 m Observations 

We targeted ah 94 EGOs in the IC08I catalog visible 
from Nobeyama (those in the northern Galactic Plane, 
6 > —20°). Our sample sources are listed in Table [U 
along with information from the literature on their MIR 
properties and CH3OII maser associations. The NII3 
(J, K)—{1,1), (2,2), and (3,3) inversion transitions and 
the 22.235 GHz H2O maser line were observed simul- 
taneously with the Nobeyama Radio Observatory 45- 
m telescope (NR045)i3 in 2008-2010. During our win- 
ter (January/February) observing sessions, the system 
temperature was typically ^100-160 K. The beamsize 
and main-beam efficiency of the NR045 at 22 GHz are 
^FWHP = 73" and ?7mb = 0.825, re spect ively. We pointed 
at the EGO positions tabulated in lC08l which are the po- 
sitions of the brightest 4.5 emission associated with 
each candidate outfiow. We note that these positions will 

The 45-m radio telescope is operated by tlie Nobeyama Radio 
Observatory, a branch of the National Astronomical Observatory 
of Japan, National Institutes of Natural Sciences. 



not necessarily be those of the driving sources (which in 
many cases a re di fficult to identify solely from the MIR 
data, see also lCOSi ). though in most cases the NRO beam 
is large enough to encompass likely driving sources as well 
as the 4.5 /im extent of the EGO. 

We used the H22 receiver, a cooled HEMT receiver, 
and eight high-resolution acousto-optic spectrometers 
(AOSs) to observe both polarizations for each line simul- 
taneously. The bandwidth and spectral resolution of the 
AOSs are 40 MHz and 37 kHz, respectively, correspond- 
ing to velocity coverage of ^500 km s~^ and resolution of 
~0.5 km s""'^ for the observed lines. The spectral chan- 
nels were Nyquist-sampled. 

The observations were conducted in position-switching 
mode, using "off" positions ~5' away. All spectra were 
checked for evidence of emission in the chosen "off" po- 
sition, and, if necessary, reobserved. Initially, each tar- 
get was observed for 2 minutes (on-source). The spectra 
were then inspected, and weak sources were reobserved 
to improve the signal-to-noise as time permitted. The 
pointing was measured and adjusted at the beginning of 
each observing run using Galactic maser sources. The 
absolute pointing of the NR045 is very accurate for 22 
GHz observations, from a few arcsec (no wind) to ^ 10" 
in the windiest conditions in which we observed-still a 
small fraction of the beamsize at 22 GHz. 

The data reduction followed standard procedu res us- 
ing t he NRO NEWSTAR software package (Iked a et al.l 
|2001|) . For each spectrum, emission- free channels were 
used to estimate and subtract a linear spectral base- 
line. For each line, the two polarizations were then 
co-added, weighted based on system temperature. The 
temperature scale was calibrated to the antenna temper- 
ature (T^) in Kelvin with the standard chopper- wheel 
method, and the main-beam temperature (Tmb) calcu- 
lated as Tmb = T^/tjub- For the H2O maser data, we 
then convert to the Jansky scale to facilitate comparisons 
with other surveys. 

Histograms of the rms are shown in Figure [1] The me- 
dian la rms is -50, 51, and 52 mK for NH3(1,1), (2,2), 
and (3,3), respectively. For our H2O maser observations, 
the median la rms is ■^O.ll Jy, corresponding to a me- 
dian 4(7 detection limit of —0.44 Jy. 



2.2. NH3 Modeling and Physical Parameter Estimation 

We estimate physical properties from the observed 
NH3 spec tra foll o wing the philosophy developed by 
Rosolowsk y et al.l (I2008D and adapted for u se in 
Dunham et al.l (|2010[ ) and IDunham et al.l (|2011b( ). The 
emission is modeled as a beam-filling slab of NH3 with 
a variable column density (A^nhs), kinetic temperature 
(Tkin), excitation temperature (Tex), Gaussian line width 
(dv), and LSR velocity (ulsr)- The model assumes the 
molecules are in thermodynamic equilibrium using an 
ortho-to-para ra tio of 1:1, which is t he high temperature 
formation limit (Takan o et al .1(200^ . Hence, the ammo- 
nia molecules are partitioned among the energy levels 
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Figure 1. Histograms of tlio distributions of rins noise for the sources in our sample for the four observed lines. 
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Here, J and K are the rotational quantum numbers 
of NH3 and, for the metastable inversion transitions, 
J = K. The energy difference, AE{J, K,i), is the split- 
ting of the symmetric and anti-symmetric states, rep- 
resenting both levels of the inversion transition. The 
antisymmetric state, AE{J,K,1), is AE/k ^ 1.1 K 
above the symmetric state {AE{J, K,0) = 0). The col- 
umn density of the molecules in the A^nh3(>/, K, i) ortho 
state is thus N]^ii3Zo{J,i)/{'2Zo) and in the para state 
N-^-RsZpiJji)/ {2Zp), where the factor of two arises be- 
cause of the assumption of a 1:1 ortho-to-para ratio. 

The optical depths in the individual transitions are 
calculated from the column densities in the individual 
states. The optical depth, hyperfine structure, veloc- 
ity information and excitation conditions are then used 
to model the individual spectra. Free parameters are 
optimized using the MPFIT least-sq uares mi nimization 
routine including parameter bounds (IMarkwa rdt 2009). 
Uncertainties in the derived parameters are also deter- 
mined from this optimization, accounting for the covari- 
ance between the parameters. We note that parameter 
uncertainties cannot account for systematic errors stem- 



ming from the uniform slab model being an incomplete 
description of the physical system. In all cases, derived 
quantities should be considered summary properties of 
the system and not a complete description. In most 
cases, this simple model reproduces the emission features 
observed on the large scales sampled. 

For some sources in our sample, however, a single-slab 
model does not adequately represent the amplitudes of 
all three NII3 transitions. Figure [5] shows examples of 
the two cases that prompted a revision of our model: (1) 
spectra that showed velocity components with different 
^LSR or (Tv; and (2) spectra that could not be well repre- 
sented by a single temperature fit. We found that includ- 
ing a second comp onent produced sign ificantly better fits 
in these cases (see iRosolowskv et al.l 2008. for more de- 
tails). A second component was introduced for any fit 
where the psr degree of freedom was larger than two 
for any individual inversion line (23 sources, ~25% of 
our sample). For two sources that met this criterion, 
the best-fit two component model included a component 
with an unphysically low excitation temperature (<2.73 
K). For these sources (G14.33-0.64 and G19.36-0.03), 
we retain the single component fits (leaving 21 sources 
with two-component fits). In the two component model, 
the two slabs are nominally beam-filling, but no radia- 
tive transfer is performed from one slab through the 
other. We see no evidence for absorption of one com- 
ponent through the other in the spectra, suggesting such 
a treatment is not needed. A simple two-component fit 
yields a substantial improvement in the quality of the 
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Figure 2. Single component (left) and two component (right) fits to sample NH3 spectra. The best-fit models are overplotted on the 
observed spectra. For the two-component fits, model spectra for each component are shown (dashed line: warmer component; dotted line: 
cooler component), as well as their sum (solid line). The 'D' at upper right in each panel indicates that our 4cr detection criterion was met 
for that transition. Gf 2.91— 0.26 (top) has two velocity components. For Gil. 92— 0.61 (bottom), two temperature components significantly 
improve the fit to the NH3(1,1), (2,2) and (3,3) spectra. 



fit for many sources, successfully identifying two veloc- 
ity/temperature components. We again note, however, 
that slab models are an incomplete description of the 
physical system; the best-fit physical parameters of the 
two components are thus likely representative but not 
definitive. We also note that a contradiction arises be- 
cause the model takes Tmb — VffiTex ~ Thg){l — e^'^) 
where r]ff — 1 is the assumed beam filling factor. How- 
ever, the parameter rjjf is degenerate with T^x, and our 
assumption that r/fj = 1 means T^^ is a lower limit. 
Relaxing this constraint owrjff leaves T^x undetermined 
for the two components, and suggests that the success 
of the simple two component fitting means the two NH3 
components are spatially distinguished on smaller scales. 

3. RESULTS 

3.1. Detection Rates 

3.1.1. Water Masers 

We define a water maser detection as >4ct emission 
in at least two adjacent channels. The overall detection 
rate is 68% (64/94), and Table H summarizes the H2O 
maser detection rates towards various EGO subsamples. 
The uncertainties quoted in Table [2] were calculated us- 



ing binomial statistics. Throughout, we treat each EGO 
separately, though we note that for EGOs separated on 
the sky by <36.5"(half the FWHP Nobeyama beam), 
our data are insufficient to determine whether one or all 
are associated with H2O masers. An unavoidable limi- 
tation of single-dish surveys is the possibility that some 
H2O maser detections are chance alignments within the 
single-dish beam, and not physically associated with the 
target EGOs. While this can only be definitively ad- 
dressed by future high-resolution observations of all de- 
tected EGOs, available data suggests the effect on the 
sample as a whole is small. We searched the literature 
for reported H2O masers with interferometric positions 
within 2' of each EGO with a H2O maser detection in our 
survey. Of 27 sources with such data available, there are 
only 3 cases ('^11%) of H2O masers within the Nobeyama 
beam and not associated with the EGO (see also §3.3p . 

One of the goals of this survey is to in vestigate whether 
the MIR EGO classifications from lCOSi correspond to dif- 
ferences in H2O maser associations or dense gas proper- 
ties. We find a somewhat higher H2O maser detection 
rate for EGOs classified as 'likely' MYSO outflow can- 
didates, compared to those classified as 'possible' based 
on their MIR properties. Two-tailed binomial tests re- 
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ject the null hypothesis that these two detection rates are 
the same at the 5% significance level (p-values ■~0.02). 
We also find a slightly higher H2O maser detection rate 
towards EGOs not associated with IRDCs, compared to 
EGOs that are associated with IRDCs. In this case, how- 
ever, two-tailed binomial tests are consistent with the de- 
tection rates being the same, at the 5% significance level 
(p-value=0.07(0.10) adopting the non-IRDC(IRDC) de- 
tection rate as the null hypothesis). If, instead, EGOs 
are grouped based on the NII3 transitions detected in our 
survey, much larger differences in the II2O maser detec- 
tion rates emerge. We detect H2O masers towards only 
44% of EGOs with NIl3(l,l) emission only, compared to 
81% of EGOs with emission in the higher-excitation NII3 
transitions: a difference of nearly a factor of two. 

There are comparably striking differences in the II2O 
maser detection rates towards EGO subsamples defined 
based on CII3OH maser associations (see Table [2]). To 
group EGOs by their CH3OH maser associations, we use 
the data in Table 1 of iCElll . This dataset, derived from 
single-dish surveys, is the most uniform available that in- 
cludes the majority (~3/4) of our northern EGO targets. 
iGEli searched for 95 GHz Class I CH3OH masers to- 
wards 192 EGOs (northern and southern) with the MO- 
PRA telescope (0fwhp~36", 3cr~1.6 Jy). They also ob- 
served EGOs without known Class II masers at 6.7 GHz 
with the University of Tasmania Mt. Pleasant telescope 
(^FWHP~7', Jy). This produced a three-tiered 

classification for Class II maser associations: (1) EGOs 
associated with Class II masers, based on published high- 
resolution data (maser positions known to ^1" or bet- 
ter); (2) EGOs for which 6.7 GHz emission was detected 
in the large Mt Pleasant beam but no positional infor- 
mation was available ("no information"); and (3) EGOs 
undetected in the Mt Pleasant observations. For this 
re ason, d efinitive Class II maser information is available 
in ICElll for a smaller number of the EGOs in our sam- 
ple (51) than for Class I masers (69 EGOs). We note 
one addition al cav eat. The 95 GHz Class I transition 
observed by ICElll is generally weaker than that at 44 
GHz, and their MOPRA observation s are significantly 
less sensitive th an th e VLA survey of [col. As a result, 
one source in the lC09l sample that has weak 44 GHz Class 
I iriasers (G37.48— 0.10) is listed as a Class I nondetection 
in CEll. 

The most dramatic difference in H2O maser detection 
rates in our survey is between EGOs associated with 
both Class I and II CH3OH masers (20/21 -95%) and 
EGOs associated with neither type of CH3OH maser 
(5/15 ~33%). The H2O maser detection rate is also very 
high (~90%) towards EGOs with Class I CH3OH masers 
(considering all Class I detections, regardless of Class II 
detections/information). This correlation is consistent 
with both Class I CH3OH and H2O masers being associ- 
ated with outflows, though H2O masers are also observed 
towards ~46% of EGOs undetected at 95 GHz. Unfortu- 
nately, comparison of "Class I only" and "Class II only" 
EGO subsamples is limited by the small number statis- 
tics. 

3.1.2. NH3 

The vast majority (97%) of our target EGOs are de- 
tected in NH3(1,1) at the 4a level (peak/rms). For a 
significant fraction (34%) of our sample, (1,1) is the only 
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Figure 3. Observed NH3 spectra with best-fit single component 
model overlaid. A 'D' in the upper right corner of a panel indicates 
that our Aa detection criterion was met for that transition. A 
complete figure set, including all 91 EGOs detected in NH3(1,1), 
is available in the online journal. 



NH3 transition detected. As shown in Tabled it is the 
detection rates for the higher-energy transitions, partic- 
ularly (3,3), that show significant differences across EGO 
subsamples. The NH3(3,3) detection rate towards EGOs 
associated with IRDCs is about twice that for non-IRDC 
EGOs; similarly, the detection r ate t owards "likely" out- 
flow candidates (as classified bv lCOSD is about twice that 
for "possible" outfiow candidates. The (2,2) detection 
rates show the same trends. 

The strongest correlation we see, however, is again with 
CH3OH maser associations. The highest (3,3) detection 
rate of any subsample is 81%, towards EGOs with both 
Class I and II masers, while the lowest (7%) is towards 
EGOs without either type. The detection rate towards 
EGOs with Class I masers (regardless of Class II associ- 
ation/information) is similarly high, at 76%. The (2,2) 
and (3,3) detection rates show similar trends, with (3,3) 
showing larger differences between subsamples. 

3.1.3. NH3 Nondetections 

Our extremely high NH3(1,1) detection rate raises the 
question of whether the three nondetections are in some 
way unusual, or interlopers in the EGO sample. The 
(1,1) nondetections do have some common characteris- 
tics: they are not associated with IRDCs and do not 
have Class I CH3OH masers. Two have detected H2O 
maser emission in our survey. G49. 424-0. 33, a IC08I 
"likely" outfiow candidate, was included in the CO!^ sam- 
ple and detected in thermal HCO+(3-2), H"CO+(3-2), 
and CH3OH (52,3-4i,3) emission with the JCMT. Thus, 
there is dense gas associated with the EGO: in combi- 
nation with the detection of Class II CH3OH (£09) and 
H2O masers (Tabled, strong evidence for the presence 
of MYSO(s). This EGO is among the most distant in 
our sample, so our Nobeyama NH3 nondetection may be 
attributable to sensitivity and/or beam dilution. 

We al so de tect H2O maser emission towards G53.92- 
0.07, a IC08I "possible" outflow candidate. Its MIR 

One source meets our 4cr peak/rms detection criterion for 
NH3(1,1) and (3,3), but not (2,2). This could be indicative of 
nonthermal (3,3) emission; however, the NH3 emission is weak and 
the (3,3) detection is marginal (<5a-). Thus, we conservatively 
treat this source as an NH3 (l,l)-only detection in our analysis of 
detection rates and in Table [2| 
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Figure 4. Histograms showing distributions of observed NH3 properties and physical properties obtained from the NH3 modeling. Bin sizes are 
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EGOs detected in NH3(1,1) are included in the first five panels ((1,1) peak, cTv, T(i^i), Tjff, and column density). Sources for which Tex—Tkin (the 
upper limit, for i^ff — l) are excluded from the filling fraction plot. EGOs detected in both NH3(1,1) and (2,2) are included in the T^^^ and (2,2) 
peak histograms, and EGOs detected in all three NH3 transitions are included in the (3,3) peak plot. 
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Figure 5. Observed NH3 spectra with best-fit two component 
model overlaid. Model spectra for each component are shown 
(dashed line: warmer component; dotted line: cooler component), 
as well as their sum (solid line). A 'D' in the upper right corner 
of a panel indicates that our 4(t detection criterion was met for 
that transition. A complete figure set, including all 21 EGOs with 
two-component fits, is available in the online journal. 

morphology is unusual amongst the EGO sample; the 
"green" source appears embedded in an 8 /im-bright pil- 
lar, and the 4.5 /xm emission is only slightly extended. 
Little is known about this source beyond its identifica- 
tion as an EGO and its association with a BGPS 1.1 
mm source, but it is possible it may be a comparatively 
evolved outlier in the EGO sample. 

Finally, G57. 61-1-0. 02 is a "possible" outflow candidate 
located on the edge of an 8 and 24/im-bright nebula, 
likely a more evolved source (e.g. compact or UC HII 
region). Formally undetected by our 4cr criteria, we do 
see a weak {^3.9a) NH3(1,1) line in our spectra (see also 

3.2. NH3 Properties 

Table [3] presents the physical properties obtained from 
the single-component NH3 modeling for all EGOs de- 
tected in NH3 emission in our survey. The NH3(1,1), 
(2,2), and (3,3) peaks (Tmb) are also listed, with 4(t 
upper limits given for undetected transitions (for all 
sources, including NH3 nondetections). If NH3(2,2) is 
not detected, the best-fit Tkin is treated as an upper 
limit and is indicated as such in Table [31 The observed 
NH3 spectra for each detected source, overlaid with the 
best-fit model, are shown in Figure [3] (available online in 
its entirety), and the property distributions for our EGO 
sample are shown in Figure HI Throughout, the NH3(1,1) 
peak (Tmb), fv, '''(1,1), ??//, and NH3 column density 
are presented for all EGOs detected in NH3(1,1) emis- 
sion. In Figure m the Tkin and NIl3(2,2) peak distribu- 
tions include only sources with >4cr NIl3(2,2) detections, 
and the NH3(3,3) peak distribution includes only sources 
with >4(T detections in all three NH3 transitions. For 
EGOs where a two component model provides a better 
fit to the observed NII3 emission f §2.2|) . Figure [5] (avail- 
able online in its entirety) shows the spectra overlaid with 
the best- fit two component model, and Table H] presents 
the parameters of the two-component fits. 

3.2.1. Kinematic Distances 

We calculate kinematic distances based on th e NH3 
velocities in Table |3] and the prescription of .Reid et al.l 
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(|2009l) . using updated input parameters (M. Reid, priv. 
comm., 2012; Galactic: Ro= 8.40 kpc, 60= 245.0 
km s-\ de/dr^ 1.0 km s^^ kpc^^; Solar: Uo= 10.00 
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km s-\ Vo= 12.00 km s-\ Wo 
pecular motions: Us= 5.00 km s^^, Vs— -6.00 km s 
Ws= 0.00 km s~^; and an assumed vlsr uncertainty of 7 
km s~^). For sources with distance ambiguities, the near 
kinematic distance is listed in Table [3J unless otherwise 
noted. The angular extent of EGOs on the sky supports 
adopting the near kinematic distance, as does the asso- 
ciation of EGOs, as a population, with IRDCs (see also 
C08i; C09). In t heir HI self-absorption study of 6.7 GHz 
CH3OH masers. [Green fc McClure-Grifhthsl (poTl have 
recently suggested assigning the far distance to masers 
associated with a few (eight) of our targets. Most of 
these assignments are 'Class B' in their scheme, reflect- 
ing uncertainty in the classification. For these sources, 
we adopt the far distance calculated from the NH3 veloc- 
ity. Maser parallax distances are adopted when available, 
as noted in Table [3} 

Three sources are undetected in NH3(1,1), and so 
present special cases for calculating kinematic distances. 
For G49.42-K0. 33, we use the Hi3CO+(3-2) velocity from 
[cm (see also S [3T3| . For G53.92-0.07, the H2O maser 
emission is very narrow (Av=1.3 km s~^), and we cal- 
culate a kinematic distance using the H2O maser peak 
velocity (Table [6]). In G57.61+0.02, we detect weak 
NH3(1,1) emission at ~3.9cr, just below our formal de- 
tection limit. The fitted vlsr of 37.4±0.1 km s~^ gives 
a kinematic distance of 4.50±1.96 kpc. For complete- 
ness, we include this source in the distance histogram 
shown in Figure [6l but not in the subsequent analysis. 
The mean(median) distance for our sample is 4.3 kpc(4.2 
kpc). 

3.2.2. Comparison of EGO subsamples 

As discussed in t i3.1.21 detection rates for the higher- 
excitation NH3 transitions differ for various EGO sub- 
samples. We consider seven pairs of EGO subsam- 
ples: (1) "likely" / "possible" outflow candidates; (2) 
sources associated/not associated with IRDCs; (3) H2O 
maser detections/nondetections in our survey; (4) Class 
I CH3OH maser detections/nondetections (regardless of 
Class II association); (5) Class II CH3OH maser detec- 
tions/nondetections (regardless of Class I association); 
(6) EGOs associated with only Class I/only Class II 
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Figure 7. NH3 property distributions for EGOs classified as "likely" and "possible" MYSO outflow candidates by ICOSI . "Likely" and 
"possible" sources are plotted as horizontally and diagonally hatched histograms, respectively. Bin sizes are the same as in Figure |4] 



CH3OH masers; and (7) EGOs associated with both 
Class I and II CII3OH masers/EGOs associated with nei- 
ther CH3OII maser type. To assess whether the NII3 
properties of these subsamples exhibit statistically sig- 
nificant differences, we ran two-sided K-S tests of eight 
parameters: the NH3 (1,1), (2,2), and (3,3) peaks (Tmb), 
^v, Vffi N(NH3), and Tkin- To maximize our 

sample size, we used the parameters from the single- 
component fits. To check for biases due to sensitiy- 
ity limits, we also ran two-sided K-S tests of distance 
and the NH3(1,1) rms for the same seven EGO sub- 
samples. Table [S] lists the subsample/parameter com- 
binations that have significantly different distributions, 
adopting a moderately conservative threshold of <0.01 
for the significance of the K-S statistic. Note that K-S 
tests involving the CII3OII maser subsamples are lim- 
ited by small sample sizes, particularly for parameters 
that require (2,2) or (3,3) detections. While we ran K-S 
tests in all cases where the subsamples being compared 
each have >4 members, we interpret the small-n results 
with caution. Statistically significant differences are seen 
most often in the NIl3(l,l) peak temperature, CTv, the 
NH3 column density, and the kinetic temperature. The 
distributions of these properties for the various subsam- 
ples are shown in Figures [7lfT2l 

The most dramatic difference is between the CTv dis- 
tributions for EGOs that are/are not detected in II2O 
maser emission in our survey (Fig. [9]). The NH3 lines are 
broader towards EGOs associated with H2O masers, with 
median CTv of 1.18 km and 0.80 km s~^ for H2O maser 



detections and nondetections, respectively (uv = ^J^™ 

). This is in agreement with previous single-dish studies 
of II2O masers in star- forming regions. In their NIl3(l,l) 
survey of 164 H2O masers (6'fwhp '-^1.4^). lAnglada et al.l 
(|1996( ) found a correlation between LH2O and the NII3 
line width; comparing their data with other NII3 sur- 
veys, they found increased NII3 linewidths towards star- 
forming regions with H2O mas ers. Both our results and 
those of lAnglada et al.l ()1996D are consistent with the 
II2O masers being excited in outflows, which also con- 
tribute to gas motions in the surrounding clump, increas- 
ing the NH3 line width. Indeed, in high-resolution Karl 
G. Jansky Very Large Arra y (VLA) observations of one 
of the EGOs in our sample, iBrogan et al. I (|2011t ) detect 
a hot (220 K), blueshifted outflow component in NH3 
emission, coincident with redshifted H2O masers. In our 
survey, EGOs with H2O masers are also generally found 
in clumps with higher NII3 column densities and higher 
kinetic temperatures than H2O maser nondetections. 

The populations of EGOs associated and not associ- 
ated with IRDCs show statistically significant differences 
in three NII3 properties: NH3(1,1) peak, a^, and the 
beam filling fraction, rjff. EGOs associated with IRDCs 
have stronger NIl3(l,l) emission (higher NH3(1,1) peak 
temperatures) and narrower NH3 linewidths (Fig.|S]). We 
note that the distance distributions for EGOs associ- 
ated/not associated with IRDCs are statistically indistin- 
guishable based on our K-S tests (K-S significance 0.21, 
median distance 4.0 and 4.3 kpc, respectively; see also 
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Figure 8. NH3 property distributions for EGOs associated/not associated with IRDCs, plotted as horizontally and diagonally hatched 
histograms, respectively. Bin sizes are the same as in Figure |4] 



gXT]). iPiUai et alj (|2006bD found that IRDCs had, on 
average, narrower NH3 linewidths than /i?^5'-selected 
high-mass protostellar objects or UC HII regions. It is 
perhaps surprising, however, that we see a difference in 
the hnewidth distributions for IRDC/non-IRDC EGOs, 
since we are specifically targeting active star-forming re- 
gions within IRDCs. The effect may be attributable 
to emission from more quiescent regions of IRDCs be- 
ing included within the Nobeyama beam (73"'^ 1.4 pc at 
a typical distance of 4 kpc) . As shown in Figure [131 
EGOs associated with IRDCs also generally have larger 
(though still small, <0.2) beam filling fractions. This is 
consistent with numerous studies that show NII3 emis- 
sion overall follows 8 fim extinctio n in IRDCs, while 
exhibiting clumpy substructure (e.g. iPillai et al.l 120065 
iDevine et aT1l20lH iRagan et al.ll2011| ). 

Interestingly, there is little evidence for statistically 
significant differences between the NII3 properties of 
"likely" and "possible" outflow candidates. The only 
properties for which the K-S significance meets our cri- 
terion are the NIl3(l,l) and (2,2) peak temperatures. 
However, their significance values are close to our cutoff 
(Table [5|), and no comparable difference is seen in the 
distributions of the physical properties (N(NH3), Tkin, 
etc.). This suggests that the difference in Tmb(1,1) and 
Tmb(2,2) might not refiect intrinsic source properties. 
We find no statistically significant difference in the dis- 
tance distributions of "likely" and "possible" EGOs. Ex- 
isting data is insufficient to evaluate other possible ef- 
fects, such as the peak 4.5 /xm positions cataloged by 



|COl(and so our pointing positions) being systematically 
further from the driving sources in "possible" EGOs. 

EGO subsamples based on CH3OH maser associations 
show notable differences in H2O maser and NH3(2,2) and 
(3,3) detection rates f §3.1|) . The K-S test analysis in- 
dicates that these CH3OH maser subsamples also have 
statistically significant differences in their NII3 proper- 
ties (Table E]). EGOs as sociate d with C lass I CII3OH 
masers (in the study of iCElll ^STT]) have brighter 
NH3(l,l) emission (e.g. greater NH3( 1,1) peak tempera- 
tures), broader NII3 linewidths, and higher NH3 column 
densities and kinetic temperatures than Class I CH3OH 
maser nondetections (Fig. [TU)) . Class II CH3OH maser 
detections/nondetections show the same trends in the 
same properties (Fig. [TT|) . EGOs associated with both 
Class I and II CH3OH masers likewise show stronger 
NH3(1,1) emission and increased NH3 linewidths and col- 
umn densities compared to EGOs associated with nei- 
ther type of CII3OII maser. Too few EGOs with neither 
CII3OII maser association are detected in NIl3(2,2) to 
run a K-S test on Tkin, but Figure [T2l shows that the ki- 
netic temperature is indeed also higher towards EGOs 
with Class I and II CH3OH masers. Of the CH3OH 
maser subsamples, the most significant difference (lowest 
K-S significance) is between the N(NH3) distributions for 
EGOs with/without Class II CII3OII masers. 

The majority of our sample of Class II CH3OH maser 
detections (21/28), and about half of our sample of 
Class I CH3OII maser detections (21/41), are comprised 
of EGOs associated with both Class I and II CH3OH 
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Figure 9. NH3 property distributions for EGOs that are/are not detected in H2O maser emission in our survey, plotted as horizontally 
and diagonally hatched histograms, respectively. Bin sizes are the same as in Figure |4] 



masers. Similarly, the majority of the Class II nondetec- 
tions (15/23) and ^Xj^ the Class I nondetections (15/28) 
are EGOs with neither type of CH3OH maser. Thus, it 
is not surprising that the Class I detection/nondetection. 
Class II detection/nondetection, and both(Class I and 
II) / neither EGO subsamples show similar patterns in 
their NH3 properties. The sample sizes of EGOs known 
to be associated with only Class I or only Class II 
CII3OII masers are small (Table [2]). Nonetheless, there 
are no indications of systematic differences in the NII3 
properties of Class I-only and Class Il-only EGOs, either 
in the K-S test results or in the plots shown in Figure [T^ 

3.3. Water Maser Properties 

For each EGO with detected H2O maser emission in 
our survey. Table [B] lists the rms, peak flux density, ve- 
locity of peak maser emission, minimum and maximum 
velocities of maser emission (>4fT, see also ? I3.1.1|) . in- 
tegrated flux density, and isotropic maser luminosity. 
Spectra are presented in Figure [14] (available online in 
its entirety), with the minimum and maximum veloci- 
ties of detected maser emission plotted as dotted lines. 
In the absence of precise positions, the extreme variabil- 
ity of H2O masers makes it very difficult to establish 
with confidence whether or not a newly observed maser 
is identifiable with one pre viously reported (as discussed 
in lBreen fc Ellingsenll20lll and references therein). The 
present study is, to our knowledge, the first systematic 
search for II2O maser emission towards EGOs. We note 
in Table[6]H2 masers detected in high- resolution stud- 



ies targeting other sa mples that fall within the polygonal 
EGO apertures from IC08I . but do not attempt to cor- 
relate our Nobeyama spectra with previ ous single-dish 
detec t ions. As in similar studies (e.g. lAnglada et al.l 
119961 : lUrguhart et al.|[201lD . we estimate the isotropic 
H2O maser luminosity, L(H2 0), as 



L{H20) 



2.30 X 10" 



S S,dV 
Jy km s^ 



D 
kpc 



(3) 



where D is the distance to the source fi )3.2.1[ Table[3l) and 
J SudV ^ is calculated over all channels that 

meet our Aa detection criterion. For H2O maser nonde- 
tections. Table [7] lists the rms and upper limit for the 
isotropic II2O maser luminosity (calculated from equa- 
tion 01 for 4cr and two channels). The distributions of 
II2O maser peak and integrated flux, luminosity, and ve- 
locity range for II2O maser detections in our sample are 
shown in Figure fTSl 

3.3.1. High Velocity Features 

H2O masers are known for their wide velocity ranges 
and high- velocity features, as compared to other masers 
found in MSFRs (e.g. CH3OH and OH). The velocity 
of the strongest II2O maser emission in a given source 
is nonetheless gen erally well-correlated w i th the vlsb. of 
the dense gas (e.g. Churchwe ll et al.lll990HAng lada et al.l 
[1996; .Urauhart et al...2011. ). Notably, for the Red MSX 
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Figur e 10. NH3 property distributions for EGOs that are/are not associated with Class I CH3OH maser emission (in ICElH . see also 
^3.1.11 1. plotted as horizontally and diagonally hatched histograms, respectively. Bin sizes are the same as in Figure l4l 



Source (RMS)!^ sample of MIR-bright MYSOs and UC 
HIIs, the distribution of Vf/aO.peafc — Vnhs is skewed to- 
wards negative velocities. The offset (from zero) is statis- 
tically significant, and indicates that blueshifted masers 
are stronger and more prevalent than redshifted masers 
(jUrguhart et al.|[2011[ ). In our sample of 62 sources de- 
tected in both H2O maser and NH3(1,1) emission, the 
mean offset VH20,peak ~ VnHs is —2.43 km s^^ and the 
median offset —0.54 km s^^. However, in our (smaller) 
sample, the offset from zero is not statistically significant 
(standard errors 1.37 and 1.72, respectively). The distri- 
bution of Vf/aCpeafe ~Vnh3 for our EGO sample is shown 
in Figure [TBI 

The relative frequency of blue- and red-shifted emis- 
sion can also be accessed by examining high veloc- 
ity maser featu res (generally defined as V — Vlsr > 
30 km s~^, e.g. iCaswell fc B rccn 2010': Urguhar t et al.l 
I2011D . ICaswell fc Breed (|201G ) recently analyzed high- 
velocity emission in numerous H2O maser subsamples 
and proposed that an excess of sources showing only 
blueshifted high-velocity emission is an indicator of 
youth. For H2O masers associated with Cla ss II CH3OH 
but no t OH masers (from the sample of iBreen et al.l 
l2010bD . they find a "blue" (blueshifted high velocity 
emission only) fraction of 16%, a "red" fraction of 8%, 
and a "red-|-blue" (both blue- and redshifted high ve- 
locity H2O maser features) fraction of 7%. Interestingly, 

For additional details o n the RMS sample, see IHoare et al.l 
ll2005l) : rUrQuhart et al.l l l200l) . 



lUrguhart et al.l ()2011t ) find a similar ratio of "blue" :"red" 
sources in their much larger sample of RMS YSOs and 
UC HII regions, though a smaller overall fraction (22%) 
of their detected H2O masers show some high velocity 
emission. Twelve of our EGO targets (~19%) have high 
velocity H2O maser features (offset by > 30 km s~^ from 
theNHg vlsr)- 6 "blue", 1 "red", and 5 "red-|-blue" . Of 
these 12 EGOs, 5 are associated with both Class I and 
II CH3OH masers, 4 are associated with Class I CH3OH 
masers and are classified as Class II "no information" 
in CEli, 1 is associated with Class I but not Cl ass II 
CH3OH masers, and 2 are not i ncluded in 'CEll'. Our 
sample sizes and those of ICasw ell & Brccn (2010) are 
too small to warrant detailed comparisons; however, the 
"blue":" red" excess we observe is generally consistent 
with their results for CH3OH maser sources. 

3.3.2. Comparison of EGO Subsamples 

To look for differences in the properties of H2O masers 
associated with the various EGO subsamples, we ran 
two-sided K-S tests for four parameters: velocity range, 
peak intensity, integrated intensity, and isotropic lumi- 
nosity. The subsample pairs considered were the same as 
outlined above ( ij3.2.2p . with the exception of H2O maser 
detections/nondetections (since we are investigating H2O 
maser properties, only detections are considered). We 
find no evidence for statistically significant differences. 
As an example, Figures ITTlfTSI show histograms of the 
H2O maser luminosity, shaded by subsample, for the six 
subsample pairs. 
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Figur e 11. NH3 property distributions for EGOs that are/are not associated with Class II CH3OH maser emission (in ICElH see also 
^3.1.11 1. plotted as horizontally and diagonally hatched histograms, respectively. Bin sizes are the same as in Figure l4l 



3.4. Properties of Associated Dust Clumps 

Of the 94 northern EGOs in our survey, 82 fah 
with the coverage of the 1.1 mm Bolocam Galactic 
Plane Survey (BGPS, re solution 33": IAguirre et al.ll2011l: 
iRosolowskv et al.l I2O1O0 . and 77 are associated with 
BGPS sources. The BGPS source extraction algo- 
rithm, Bolocat, uses a seeded watershed approach to 
identify the boundaries of BGPS sources, and outputs 
'label maps' in which each pixel assigned to a source 
has a val ue of that sou r ce's BGPS c atalog number (see 
[Rosolowskv et alj I201CH : IDunham et al.i i2011a . for more 
details) . If the position of an EGO from lCOSl falls within 
the Bolocat-defined boundary of a BGPS source, we con- 
sider the EGO and BGPS source to be associated. 

We calculate clump gas masses from the 1.1 mm dust 
continuum emission 



4.79 X 10-^^ RS^iJy)D^{kpc) 

B(v,Tdust)Hv 



(4) 



where is the integrated flux density from the BGPS 
catalog corrected by the recommended fac tor of 1.5±0.15 
(|Aguirre et al.ll201lHD{^nham et al.ll2010D , D is the dis- 
tance to the source f ij3.2.11 Table [3]), ^{v.Tdust) is the 
Planck function, R is the gas-to-dust mass ratio (as- 
sumed to be 100), and is the dust mass opacity co- 



T he slight difference from the statistics in IDunham et al.l 
H2011al ') is because we consider G19.01— 0.03 as a single EGO, while 
they treat this EGO and its northern and southern outflow lobes 
(for which separate photometry is given in , COS,') as three objects. 



efficient in units of cm^ g~^. We follow recent BGPS 
studies (e.g. IDunham et al.l 120101 l2011airbl ) in adopting 
'«27iGHz/R.=0.0114 cm^ g~^. Our NH3 observations pro- 
vide a measurement of the clump-scale gas kinetic tem- 
perature, Tfei„, and we assume T^dust='^kin in calculat- 
ing the clump masses. To estimate the volume-averaged 
number densities of the clumps, we use the clump gas 
mass from equation|3]and the de convolved angular sourc e 
radius from the BGPS catalog (jRosolowskv et al]|2010D . 
assuming spherical geometry. For consistency with 
iHill et al] (2005) (sec ij4.2[) . we adopt a mean mass per 
particle /i =2.29 m^f. The 1.1 mm flux densities, radii, 
gas masses, and volume-averaged number densities for 
the clumps associated with our target EGOs are listed in 
Table [5] For the three BGPS sources in our sample that 
could not be stably deconvolved (listed as "null" radii in 
the BGPS catalog), we adopt half the BGPS beamsize 
as an upper limit to the source radius, e.g. R<16'.'5. The 
derived number densities for these sources are thus lower 
limits, and are indicated as such in the tables and fig- 
ures. We regard this radius upper limit as conservative 
because source radii can sometimes be determined for 
source diameters smaller than a beam width. However, 
given the substantial uncertainty in relating an emission 
distribution to a true radius, particularly at low signal- 
to-noise, a more aggress ive limit could be incorrect (e.g. 
IRosolowskv "eFaLl 120101: IRosolowskv fc Lerovl [20061 ) . If 
we instead adopted an upper limit of half the BGPS 
beamsize for the source diameter, this would increase 
the density limits by a factor of 8. 
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Figure 12. NH3 property distributions for EGOs associated with both Class I and II CH3OH masers (green), only Class I CH3OH masers 
blue) , only Class II CH3OH masers (red), and neither type of CH3OH maser (orange) (CH3OH maser associations from [CTj^ 111 , see also 
13.1.11 1. Bin sizes are the same as in Figure |4] 
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Figure 13. Distribution of the beam filling fraction, for 
EGOs associated/not associated with IRDCs, plotted as horizon- 
tally and diagonally hatched histograms, respectively. The bin size 
is 0.01. Sources for which Tex ='7kin a-iid Vff =1 ^i'g not shown. 



To estimate clump parameters consistently for the 
largest possible number of sources in our sample, we first 
calculate Mgas and n//2 as described above using the gas 
kinetic temperatures derived from the single-component 
NH3 fitting. For EGOs undetected in NH3(2,2), we treat 
the best-fit Tfcj„ as an upper limit (see also ij3.2p : the 
derived clump mass and density are thus lower limits. 
The clump masses estimated using well-determined ki- 
netic temperatures are in the range of hundreds to thou- 
sands of solar masses (Fig. [19]), with a mean (median) 



Figure 14. H2O maser spectrum. The minimum and maximum 
velocities of detected maser emission (> 4(T, Table (6)1 are shown 
as dashed vertical lines. The velocity range shown for each EGO 

extends from Ymin.water-SO km S"^ to Vmaa;,mater+30 km S"^. 

A complete figure set including spectra for all EGOs detected in 
H2O maser emission is available in the online journal. 

of -1850 M0 (~1010 Mq). The range of EGO dust 
clump masses is consistent with expectations for MYSOs 
bas ed on bolometer st udies of other samples. For exam- 
ple, [Rat^orneler^L] ([2006) find a median IRDC mass of 
-940 M0 (range -120 to 16,000 Mq), and lMueller eTall 
(|2002f ) report a similar range and a mean mass of 2020 
Mq for a sample of H2O maser sources with high lumi- 
nosities (Lboi >10^Lq). The star- forming sources (those 
with Class II CII3OH masers and/or UC HII regions) in 
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Figure 15. Histograms showing the distributions of H2O maser properties for all H2O maser detections in our sample. The panels show 
(clockwise from upper left): peak flux density, integrated flux density, velocity range of maser emission (>4o'), and isotropic H2O maser 
luminosity (S I3.3|l . Bin sizes are 0.5 dex (peak and integrated fluxes and luminosity) and 10 km s~^ (velocity range). 



the iHill et all (|2005 f) dust clum p sample similar ly span 
a mass range of ^lO^-lO'^ Mq (jHill et aUlMol ). Only 
one EGO appears to be a potential example of a nearby, 
low-mass YSO based on the properties of its associated 
dust clump: G49.91+0.37, which has a low (<10 M©) 
lower-limit mass and a near kinematic distance of 0.53 

-0.53 '^P'-- 

The significantly improved fits obtained with two tem- 
perature components for ^1/4 of our NH3 spectra in- 
dicate emission from both warmer inner regions and 
cooler outer envelopes along our lines of sight. As 
noted in H2.21 the beam filling factor r/ff and the ex- 
citation temperature T^x are degenerate for the two- 
component modeling. To estimate the relative contri- 
butions of the warm and cool components, we assume 
Tkin—Tpx and calculate rjff = ^ex-z^-TS^ each com- 

and 



-+i 



ponent. We then assign weights, Wwarm = - 
M^cooi = '7warL i i , sud rccalculatc the clump mass as 

Icool 

Mtotal = Mgas,warm + Mgas,cool whcrC 



M. 



4.79 X 10-^'^RS,{Jy)W^arn.D%kpc) 



gas.warvrL 



and 



.cool 



4.79 X lQ-^^RS,{Jy)W,ooiD\kpc) 



dust, I 



ol)^u 



(5) 



(6) 



The volume-averaged number density is then estimated 
as described above, using Mtotai in place of the single- 
temperature isothermal gas mass calculated from equa- 
tion m We can estimate revised clump masses and 
number densities in this way for 16 of the 21 sources 
with two-component NH3 fits. For these sources, the 
median(mean) mass fraction in the warm component is 
5.5%(10.0%). For five sources with two-component fits, 
Tex— T kin for one of the modeled temperature compo- 
nents (the upper limit). Coincidentally, three of these 
five sources fall outside the BGPS survey area. For the 
remaining two sources, we retain the isothermal masses 
and densities in our analysis. 

3.4.1. BGPS 1.1 mm Nondetections 

Young, actively accreting MYSOs are expected to be 
still embedded in their natal clumps; as discussed above, 
we find a strong correlation between EGOs and BGPS 
1.1 mm dust sources. The five EGOs within the BGPS 
survey area but not matched to a BGPS source are all 
detected in NH3(1,1) emission; as a group, they are not 
particularly distant (all have D< 5.5 kpc. Table [3]). The 
NH3(1,1) detections indicate that dense gas is present; 
here we briefly consider the nature of these EGOs and 
the reasons for their lack of counterparts in the BGPS 
catalog. 

The rms noise of the BGPS survey varies with Galac- 
tic longitude, a nd is locally increased in the vicinity of 
bright sources (jAguirre et al.. ,2011.) . Two of the un- 
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Figure 16. Distribution of the difference between the NH3 vi^sR 
and the velocity of peak H2O maser emission for all sources with 
both H2O maser and NH3(1,1) detections. The bin size is 2 km s^^ . 

matched EGOs (G62. 70-0.51 and G58.09-0.34) are at 
/ 60°, where the noise is signific antly higher than i n 
most of the inner Galaxy (Fig. 11 of lAguirre et"aII[20TTh : 
G62. 70-0.51 is also near the edge of the BGPS map. A 
third unmatched EGO, G49.27-0.32, is in a region of 
locally high noise due to its proximity to W51. Thus, 
it is possible that these EGOs are associated with mm 
dust clumps that would have been detected elsewhere in 
the BGPS survey. The distances of G49. 27-0.32 and 
G62.70— 0.51 are typical of our sample (D=5.5 and 3.9 
kpc respectively, Table [3|), and their properties are gen- 
erally consistent with EGOs detected only in NH3(1,1) 
emission and matched to BGPS sources. The increased 
noise of the BGPS survey at the locations of these EGOs 
thus seems to be a likely explanation for their lack of 
BGPS counterparts. G58.09— 0.34, however, may be an 
example of a nearby, low-mass YSO: it has a near kine- 
matic distance of 0.74to gi kpc and exceptionally narrow 
NH3(1,1) emission (crv-0.23). 

Examining the BGPS images suggests that the two 
other unmatched EGOs (G50.36-0.42 and G29.89-0.77) 
are associated with 1.1 mm emission, despite not be- 
ing matched to BGPS catalog sources. G50.36— 0.42 ap- 
pears to be associated with faint 1.1 mm emission that 
fell below the threshold for extra ction as a BGPS source 
(|Rosolowskv et al.ll2010[ ). A IC08I "possible" outflow can- 
didate (D=3.0 kpc), G50.36-0.42 also has detected H2O 
maser emission in our survey. G29.89— 0. 77 is immedi- 
ately adjacent to a BGPS source, but the IC08I position 
falls outside the BG PS so urce boundary defined by the 
label maps. Also a IC08I "possible" outfiow candidate, 
G29.89— 0.77 has the strongest NH3 emission of the un- 
matched EGOs; though the (2,2) line is formally unde- 
tected by our 4cr criteria, weak NH3(2,2) emission is evi- 
dent in the spectrum. Taken together, this evidence sug- 
gests G29.89-0.77 and G50.36-0.42 are likely similar in 
nature to EGOs that are matched to BGPS sources. 

4. DISCUSSION 

4.1. EGOs in Context 

4.1.1. Comparison with Other Samples 

A notable feature of EGOs, compared to other sam- 
ples of young massive (proto)stars, is their very strong 
association with both Class I and II CH3OH masers, re- 
flected in notably high detection rates in CII3OII maser 



surveys to date (e.g. [COl : ICEll . Since II2O maser and 
NH3 observations are common tools for studying mas- 
sive star formation, our Nobeyama survey allows us to 
better place EGOs in their broader context, by com- 
paring their molecular environments to those of MYSOs 
selected using other criteria/tracers. Table |9] summa- 
rizes H2O maser and NIl3(l,l) detection rates towards 
a variety of MYSO samples from the literature, chosen 
to cover a range of sample selection criteria, survey pa- 
rameters, and proposed evolutionary state of the target 
objects. The strong correlation of EGOs with 6.7 GHz 
CII3OII masers and dust clumps ( §3.4p suggests these 
as natural comparison s amples (indeed, the sam ples of 
iBreen fc Eningsen|[20lTI: iBartkiewicz et al.ll201ll include 
so me EGOs, see al s o disc ussion therein). 'Active' cores 
in [Chambers et all (j2009() are defined by the presence of 
"green fuzzy" and 24 /im emission. They define "green 
fuzzy" broadly, compared to C08 EGOs; still, one might 
expect these sources to be similar to EGOs associated 
with IRDCs. In contrast, MYSO and UC HII samples 
compiled using the IRAS or MSX point source catalogs 
comprise sources that are more MIR-bright than EGOs, 
and so li kely more luminous and/or more evolved (see 
also [cosh. 

As illustrated by Table [SI H2O maser detection rates 
towards massive (proto)star samples span a broad range, 
from <20% to >80%: our overall detection rate of 68% 
is towards the upper end of this range. Notably, our 
II2O maser detection rate towards EGOs associated with 
both Class I and II CII3OII masers (95%) exceeds, to 
our knowledge, any reported in the literature. Our much 
lower detection rate towards EGOs with neither CII3OH 
maser type (33%) is nonetheless higher than towards 
quiescent dust clumps or IRDC cores. In general, the 
II2O maser associations of EGO subsamples are simi- 
lar to those of the most comparable subsamples in Ta- 
ble [HI For example, our detection rate for EGOs associ- 
ated with Class II masers (regardless of Class I associa- 
tion) is roughly comparable to those for Class II CII3OII 
maser and dust clump/Class II CH3OH maser samples. 
Sensitivity is of course an important consideration, par- 
ticularly in light of recent evidence that II2O maser fiux 
density increases as sources e volve, then turns over a t 
a late (UC HII region) stage ()Breen fc Ellingsenl[20Tll) . 
While H2O masers are variable, the fact that we fail to 
detect H2O maser emission towards EGO Gll.ll— 0.11, 
where a weak (^0 .3 Jy) H2O maser was reported by 
[Pillai et al.l (|2006a[ ). indicates that some EGOs are as- 
sociated with H2O masers below the detection limit of 
our survey. Most of the surveys in Table [S] have sensitiv- 
ity comparable to or better than our Nobeyama data. 

The properties of the H2O masers detected towards 
EGOs are typical of H2O masers detected towards 
MYSOs. For example, the distributions of the veloc- 
ity range of detected masers and of the velocity offset 
between dense gas and peak maser emission (Fig. [TMTOl 
see also i j3.3.1[) are generally similar to those reported 
in the literature, including for more evol ved UC HII re- 
gion s amples (e.g. |C hurchwc ll ct al. 1990; lAnglada et al.l 
119961 : [Urquhart et al.l 120111 ). Based on their study of 
MIR-bri ght MYSOs and UC H II regions from the RMS 
sample, [Urquhart et al.l (|2011[ ) argue that H2O maser 
properties (in particular, Ljso) are driven by the bolo- 
metric luminosity of the central MYSO (see also ^A.2\ . 
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Figure 17. Distributions of tiie H2O maser luminosity for different EGO subsamples. Upper left: Divided by association with IRDCs. 
EGOs associated and not associated with IRDGs are plotted as horizontally and diagonally hatched histograms, respectively. Upper right: 
Divided by "likely" /"possible" outflow candidates. EGOs classified as "likely" and "possible" by COS are plotted as horizontally and 
diagonally hatched histograms, respectively. Lower left: Divided by Class I ClfsOIf maser association (regardless of Class II associa- 
tion/information). Class I detections/nondetections are plotted as horizontally and diagonally hatched histograms, respectively. Bottom 
right: Divided by Class II CH3OH maser association (regardless of Class I association/information). Class II detections/nondetections 
are plotted as horizontally and diagonally hatched histograms, respectively. The bin size is 0.5 dex, as in Figure [Tsl The significance of 
the K-S statistics (low values indicate different cumulative distribution functions) are 0.98 (IRDC/no IRDC), 0.06 (likely/possible), 0.51 
(Class I/no Class I), and 0.51 (Class Il/no Class II), indicating no statistically significant differences in the distributions of the H2O maser 
luminosities. 
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Figure 18. Distribution of H2O maser luminosity for EGOs asso- 
ciated with both Class I and II CH3OH masers (green), only Class 
I CH3OH masers (blue), only Class II CH3OH masers (red), and 
neither type of CH3OH maser (orange). The bin size is 0.5 dex. 



The distributions of H2O maser peak and integrated flux 
density and isotropic luminosity for the MIR-bright RMS 
samp le have high-end tails (e.g. Fig. 8 of lUrguhart et ahl 
1201 If) : the strongest RMS water masers are several or- 
ders of magnitude brighter and more luminous than the 
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Figure 19. Distribution of clump masses estimated from 1.1 mm 
dust continuum em ission for sources with well-determined kinetic 
temperatures f i|3.4ll . For clarity, only nominal mass values from 
Table[8]are plotted. The bin size is 0.2 dex. 



strongest water masers we detect towards EGOs. How- 
ever, two-sided K-S tests on these parameters indicate 
that the differences are not statistically significant (K-S 
significance 0.055, 0.249, and 0.027 for Speak, Smt, and 
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Liso, respectively). The K-S tests are consistent with the 
RMS and EGO water masers being drawn from the same 
parent distribution. 

As discussed in §3.2.21 and §3.3.21 we find evidence for 
statistically significant differences among EGO subsam- 
ples in NH3 but not in H2O maser properties. Other NH3 
studies of large MYSO samples similarly find significant 
internal variations. The mean kinetic temperature, NH3 
linewidth, and NH3 column density of BGPS sources in- 
crease with the number of associated MIR sources (al- 
beit with c onsiderable scatter, p articularly in Tkim e.g. 
Fig. 23 of IDuiiham et all 1201 llj). In the RMS sample, 
lUrquhart et al.l (|2011[ ) find that the mean kinetic tem- 
perature, NH3 column density, and NH3 linewidth are 
higher for UC HII regions than for MYSOs. Overall, the 
clump-scale NII3 properties of EGOs are roughly com- 
parable to those of oth er MYSO sampl e s. C omparing 
Figure m to Figure 4 of lUrquhart et al.l (|201l !). for ex- 
ample, the linewidth, Tkm, and N(NH3) distributions 
are broadly similar (accounting for the conversion be- 
tween (Tv and FWHM linewidth), though our sample is 
considerably smaller. The distribution of NH3 column 
density extends to lower values for EGOs than for the 
RMS sample; however, this is a beam-averaged quantity, 
and the Nobeyama beam (~73") is considerably larger 
than that of the GBT. For BGPS sources, the low end of 
the NII3 column density range (also based on GBT ob- 
servations) extends to ~1.7xl0^^, more comparable to 
our EGO results. The EGO Tkin distribution (from the 
single-component fitting, for consistency with other stud- 
ies) lacks the high temperature (>40 K) t ail seen in RMS, 
UC HII region, and even BGPS samples (lUrquhart et al.l 



IMltlDunham et al.ll2fribl : IChurchwen et al.lll99nD . The 

mean Tkin for the EGO sampl e as a whole (2 3.6 K) is 
higher than that of the iDunham et a l.l (|2011b) sample 
(17.4 K, for their 'T^ subsample' consisting of (2,2) de- 
tections) and similar to that of the RMS sample as a 
whole (-22 K). 

These general comparisons illustrate that the II2O 
maser and clump-scale NH3 properties of EGOs are con- 
sistent with their being a population of young MYSOs. 
However, we emphasize that the differences within sam- 
ples (EGOs, RMS sources, BGPS sources) are often 
as great or greater than the differences between them. 
These mira-sample differences emphasize the importance 
of studying multiple star formation tracers across wave- 
length regimes. 

4.1.2. Comparison with Star Formation Criteria 

By combining our Nobeyama NH3 data with the 
BGPS, we can also consider the dust clumps associ- 
ated with EGOs in the context of proposed star for- 
mati on thresholds. Unlike purely mm-selected samples 
(e.g. IDunham et al.l l2011b( ) , all of the clumps we con- 
sider are associated with EGOs, and thus demonstra- 
bly star-forming (many are also associated with other 
MIR sources) . Figure [20] shows a mass-radius plot for 
clumps with well-determined (non-limit) Tkin and radius, 
with the clump mass estimated assuming Tdust=Tkin 
from the single-component NH3 fits. The errors bars 
shown in Figure [201 indicate the range in radius asso- 
ciated with the distance uncertainty from Table [H and 
the range in mass associated with the combined uncer- 
tainties in the BGPS integrated fiux density, the BGPS 
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Figure 20. Clump mass vs. radius for BGPS sources associated 
with EGOs. Clump masses are calculated assuming Tduat='^kin 
from single-component NH3 fitting. Open squares indicate nom- 
inal values from Table |8] The error bars indicate the range 
in radius associated with the uncertainty in distance (Table [Sjl 
and the range in mass associated with the combined uncertain- 
ties in the I3GPS integrated flux density, the BGPS flux cor- 
rection facto r, and the distance. The st ar formation thresh- 
olds of 'Krumholz & McKcc (2008'), 'Hcidcrman et al. (2010) and 
Lada ct al. (2010), and Kauffmann & Pillai (2010) are i ndic ated as 
dot-dashed, dashed, and dotted lines respectively (see 13.411 . Only 
sources for which the T^in and radius are well-determined (non- 
limit) are plotted. H2O maser detections arc plotted in green, and 
H2O maser nondetections in red. 



flux correction factor, and the distance (see also 
The error bars do not include systematic uncertainty 
in the radius estimate due to different geometries (see 
also [Rosolowsky et al. 2010). Three proposed star for- 
mation thresho lds a r e indi cated on Figure [5D1 (1) the 
IKrumholz fc McKed ([20081 ) threshold for massive star 
formation of 1 g cm~^ ( 4788 M fD pc ^); (2) the av- 
erage of the ILada et all (|201(]| ) and IHeiderman et al.l 
(j201(l) thresholds for "efficient " star forma t ion ( 122.5 
Mp, pc~^); and (3) the "kauffmalll l fc Pillai! ([2OIOI ) and 
I Kauffmann et all ([2010D threshold for massive star for- 
mation. We refer to these as the KM, HL, and KP 
thresholds, respectively. As in the recent BGPS study 
of Dunham et al, (201lB), we scale the KP criterion of 
M(r)> 870 Mq (r/pc)i-33 to M(r)> 580 Mq (r/pc)i-33 
to account for the difference between our assurn e d dus t 
opacity and that adopted by [Kauffmann &: Pillail ([2010D . 
Adopting the nominal clump mass and radius values from 
Table [3 70% (35/50) of the sources shown in Figure [201 
exceed the KP threshold, and 76% (38/ 5 0) exc eed the 
HL threshold; as in the IDunham et al.l ([2011bl ) study 
of BGPS sources, none of our EGO clumps meet the 
KM criterion. We emphasize that the points in Fig- 
ure [20I represent average surface densities over entire 
BGPS sources, and that the BGPS and Nobeyama ob- 
servations probe large scales. At a typical distance of 
4 kpc, the 33" BGPS beam is ~0.64 pc, and the 73" 
Nobeyama beam ~1.4 pc. Interferometric observations 
of EGOs, and of other MYSOs, provide ample evidence 
for substructure (e.g. cores and (proto)clusters) and vari- 
ations in gas temperature on much smaller scales (e.g. 
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Figure 21. NH3 spectra of the EGO G35.03+0.35, in which 
IBrogan et aTl ((20lJ) detect an NH3(3,3) maser. 



iCyganowski et al.ll201 la! IBrogan et al. II2011D . 

Having placed clumps on the mass-radius plot using 
(primarily) the BGPS data, we use our Nobeyama survey 
data to look for differences in the properties of clumps 
above/below the HL and KP thresholds. As in our com- 
parison of EGO subsamples (i i3.2.2|) . we ran two-sided 
K-S tests on eight NH3 parameters (the NH3 (1,1), (2,2), 
and (3,3) peaks (Tmb), CTv, T(i i), 77//, N(NH3), and 
Tfcin). We find statistically significant differences only 
for the NH3(1,1) and (2,2) peak temperatures and the 
filling fraction f?//^^, with clumps below the HL and 
KP thresholds having lower values of these parameters. 
The K-S tests indicate no statistically significant differ- 
ences in the distributions of the physical properties Uv, 
Tkin, and N(NH3) for clumps above/below the thresh- 
olds. Interestingly, and perhaps counterintuitively, the 
H2O maser detection rates are higher for EGOs asso- 
ciated with clumps below the HL and KP thresholds 
(Fig. [201). The H2O maser detection rate is 0.74(±0.07) 
for sources that meet the KP criterion, and 0.93(±0.06) 
for sources that do not (uncertainties in detection rates 
calculated using binomial statistics). Similarly, the H2O 
maser detection rates are 0.76(±0.07) and 0.92(±0.08) 
for sources that do/do not meet the HL criterion, re- 
spectively. 

The nature of the EGOs associated with clumps that 
fall below the KP threshold requires further investiga- 
tion. The higher H2O maser detection rate towards 
clumps below the KP threshold is surprising, and the 
lack of difference in NH3 properties suggests a contin- 
uum, rather than a sharp distinction. Additionally, one 
source that fahs below the KP threshold, G24. 94-^0. 07, 
is associated with 6.7 GHz Class II CH3OH maser and 
cm continuum emission (C09,; CI lb ), both indicative of 
the presence of an MYSO. We note that the placement of 
clumps on a mass-radius plot is sensitive to assumptions 
about clump temperature structure (or lack thereof). For 
EGOs in our study fit with warm and cool components, 
the warm component constitutes a small fraction of the 
clump mass; the bulk of the material generally has tem- 
perature Tcooi single comp., and SO the isothermal as- 
sumption underestimates the clump mass ( ^3.41 Tablc|5]). 
Interferometric NH3 observations show significant tem- 
perature structure on scales within the Nobeyama beam 

We note rjff is mildly degenerate with Tmb(1;1)- 
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Figure 22. Top: isotropic H2O maser luminosity vs. volume- 
averaged number density estimated using Tfci„ from single- 
component NH3 fitting. * indicates EGOs with H2O maser de- 
tections in our survey and well-determined density estimates (e.g. 
neither Tj,j„ nor R is a limit). Filled downward-pointing triangles 
indicate Aa L(H20) upper limits for EGOs with well-determined 
density estimates that are H2O maser nondetections in our sur- 
vey. EGOs for which the estimated density is a lower limit are 
represented as open triangles: open right-facing triangles indicate 
H2O maser detections, and open downward-pointing triangles 4(t 
L(H2 0) upper limits for H2O maser nondetections. Bottom: Same 
as top, except the density estimate accounts for warm and cool 
com ponen ts when the NH3 spectrum is fit with two-components 
(see ^3A\i . 



for G35.03+0.35 (Fig. 3 of IBrogan et"an[20TTf ) . a source 
that did not require two temperature components to 
fit its Nobeyama NH3 spectrum (Fig. I^Tj) . On larger 
scales, many of the BGPS sources associated with EGOs 
(and plotted in Fig. [20l) extend beyond the Nobeyama 
beam. If isothermal clump masses for EGOs tended to be 
underestimates-due to temperature structure on small or 
large scales-this would move points up in Fig.l20| and in- 
crease the proportion of sources above the KP threshold. 
Additional data-such as NH3 maps with sufficient reso- 
lution to probe the temperature structure of the BGPS 
clumps-are needed to address this issue. Interferomet- 
ric (sub)mm observations, to resolve the dust continuum 
emission and detect individual cores, and improved con- 
straints on bolometric luminosity (e.g. from HiGal) will 
also help to clarify the nature of the driving sources. 

4.2. Correlations between H2O Maser and Clump 
Properties ? 

Over the past decades, numerous authors have in- 
vestigated possible correlations amongst clump, H2O 
maser, an d driv ing sources properties i n MYSO sam- 
ples (e.g. iChurchwell et al., .1990; .Anglada et al.l 119961: 
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Figure 23. Top: isotropic H2O maser luminosity vs. Tkin from 
single-component NH3 fitting. * indicates EGOs with H2O maser 
and NH3(2,2) detections in our survey (e.g. Tkin well-determined). 
Filled downward-pointing triangles indicate 4cr L(H20) upper lim- 
its for EGOs undetected in H2O maser emission but detected in 
NH3(2,2). EGOs undetected in NH3{2,2)-for which the best-fit 
Tkin is treated as an upper limit-are represented as open trian- 
gles: open left-facing triangles indicate H2O maser detections, and 
open downward-pointing triangles Aa L{H20) upper limits for H2O 
maser nondetections. IJottom: Same as top, except for sources 
fit with two NH3 components, Tj;j„{cool) is plotted in blue and 
Tfei„(warm) in red. 

IBreen fc Ellingsed [Mil: lUrquhart et"all [20TTI ). Re- 
cently, two studies have reported correlations between 
H2O maser luminosity and the properties of the driving 
source or surrounding clump. For their sample of ^300 
RMS sources with H2O maser detections. lUrquhart et alj 
(|20Tll) find that H2O maser luminosity is positively cor- 
related with bolometric lu minosity for both MYSOs and 
HII regions. In contrast, IBreen fc EllingsenI ()2011l ) re- 
port an anticorrelation between clump H2 number den- 
sity and H2O maser luminosity, which they attribute to 
an evolutionary effect: more evolved sources have more 
luminous water masers and are associated with lower- 
density clumps. All of these studies have combined H2O 
maser and either NH3 or (sub)mm dust continuum data. 
IBreen fc EllingsenI (|201lD . in particular, caution that the 
clump densities used in their study ffrom lHill et al.ll2005l) 
were calculated assuming a single temperature for all 
clumps, and that temperature differences could create 
the apparent density trend. Our NH3 and H2O maser 
survey, in combination with the BGPS, provides the nec- 
essary data to fully explore correlations between maser 
and clump properties, and test evolutionary interpreta- 
tions. 

Figure [22] shows that when clump densities are cal- 



Figure 24. Top: isotropic H2O maser luminosity vs. clump mass, 
assuming Tdust=Tkin from single-component NH3 fitting. * in- 
dicates EGOs with H2O maser and NH3(2,2) detections in our 
survey (e.g. Tkin well-determined). Filled downward-pointing tri- 
angles indicate 4cr L(H20) upper limits for EGOs undetected in 
H2O maser emission but detected in NH3(2,2). EGOs undetected 
in NH3(2,2)-for which the best- fit T kin is an upper limit and the 
clump mass thus a lower limit ( i|3.4l| -are represented as open tri- 
angles: open right-facing triangles indicate H2O maser detections, 
and open downward-pointing triangles 4<j L(H2 0) upper limits for 
H2O maser nondetections. Bottom: Same as top, except the mass 
estimate accounts for warm and cool components for sources fit 
with two NH3 components ( ^13.411 . 

culated for our sample using measured clump tempera- 
tures, there is no correlation between H2O maser lumi- 
nosity and clump density: the log-log plot of L(H20) vs. 
number density is a scatter plot. This remains the case 
even when accounting for the contributions of warm and 
cool gas for sources that require two-component NH3 fits 
(Fig. [22I bottom; §3.4p . The partial correlation coeffi- 
cients, computed with the distance squared as an inde- 
pendent parameter, are 0.04 and 0.06 for the one- and 
two-temperature component density estimates, respec- 
tively (only sources with H2O maser detections and non- 
limit densities are included in the calculation). These 
low values confirm that H2O maser luminosity and clump 
number density are uncorrelated in our data. 

In contrast, H2O maser luminosity is weakly correlated 
with clump temperature, as shown in Figure 1231 For 
EGOs detected in both H2O maser and NIl3(2,2) emis- 
sion, the partial correlation coefficient is 0.36 for Tkin 
derived from the single-component fits (again computed 
with the distance squared as an independent parame- 
ter). Interestingly, if we recompute the partial correla- 
tion coefficient using the Tkin of the warm component for 
sources that require two-component fits (and the single- 
component Tkin for all other sources), the value is re- 
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duced to 0.22. This is somewhat surprising, since the 
warm component traces gas nearer to, and heated by, 
the central MYSO. 

We also find a weak positive correlation between H2O 
maser luminosity and clump mass (Fig. I24|) . Calculat- 
ing clump masses assuming Tdust=T^kin from the single- 
component NH3 fits, the partial correlation coefficient 
is 0.44 (for EGOs detected in both H2O maser and 
NH3(2,2) emission, so that Tkin is well-determined). 
The calculated partial correlation coefficient is very sim- 
ilar (0.43) if the presence of two temperature compo- 
nents is accounted for when estimating the clump mass 
(> j3.4[) . A K-S test indicates no statistically signifi- 
cant difference between the mass distributions of clumps 
wit h /without H2O masers, in contr ast to earlier stud- 
ies (jChambers et a l. 2009; Br^eenjLEllingsen 2011). The 
K-S statistic is 0.26 using the isothermal clump masses 
(for EGOs with (2,2) detections and so well-determined 
Tkin, as above), and increases to 0.45 if clump masses 
are estimated accounting for the two temperature com- 
ponents. Both previous studies a ssumed dust tempera- 
tures, and [Chambers et al.l (|2009 l) found that the proba- 
bility their cores with/without H2O masers were drawn 
from the same distribution increased dramatically (by a 
factor of >50, to 0.11) if they assumed a higher tem- 
perature for active cores (compared to assuming a single 
temperature for all cores). The 6'fwhm of the BGPS data 
(~33^') is larger than th at of th e SIMBA data used by 
iBreen & Ellingsen' I'M! f^24": iHiU ct al. 2005) or the 
IRAM 30-m data used b v lChambers"et a l. (2009) (--11"; 
iRathborne et al. Il2006[) . Additi onal data (such as tem - 
per ature measurement s for t he [Chambers et al.l ()2009f ) 
and lBreen fc EllingsenI (|201lD sources) would be required 
to assess whether this difference in scale contributes to 
the difference in findings. 

Our results are consistent with the positive correla- 
tion betwe en H2O maser and bo lometric luminosity re- 
ported by lUrauhart et~all (|2011[ ) for RMS sources. In 
this picture, the key factor is the bolometric luminos- 
ity of the driving MYSO, with more luminous MYSOs 
exciting more luminous H2O masers. The observed cor- 
relations of H2O maser and clump properties (temper- 
ature and mass) are then understood in terms of the 
relationship between a clump and the massive star(s) it 
forms. The final mass of an actively accreting MYSO 
is limited by the available mass reservoir, and studies 
of more evolved sources (UC HIIs) indicate that higher- 
mass clum ps form higher-mass ( and thus more luminous) 
stars (e.g. [Johnston et al.|[2009l) . The more luminous an 
MYSO, the more energy it will impart to its environs, 
and the more it will heat the gas and dust of the sur- 
rounding clump. 

4.3. NH3(3,3) Masers 

While NH3(3,3) maser emission in a MSFR was 
first reported sev e ral d ecades ago (DR21(0H); 
iMangum fc WoottenI |1994[) . the number of known 
examples-all detected with the VLA-has remained 
smah (e.g. W51, NGC6334I, IRAS 20 126-^4106 
G5.8 9 -0.39: IZ hang & Ho 1995; Kraemer AJacksonI 
119951 IZhang et al.. ,199a: .Hunter et al.. ,2008. ). Two re- 
cent, large-scale single dish surveys each report a single 
NH3(3,3) maser candidate: a blind survey of 1 00 deg^ 
of the Galactic plane (HOPS: IWalsh et alJl20Tll) . and a 



targe ted survey of ^^600 RMS sources ( Urguhart et al 



mThis paucity of candidates led lUrquhart et al 
to suggest that bright NH3(3,3) masers are rare. 
One of our targets, G35. 03-1-0. 35, wa s recently ob- 
serve d in NH3(l,l)-(6,6) with the VLA ( Brogan et all 
|2011| ). In addition to complex thermal NH3 emission 
from a (proto)cluster, nont hermal NH3(3,3) an d (6,6) 
emission are clearly detected ()Brogan et al. II20T1I Fig. 2; 
peak (3,3) intensity <70 mJy beam~^). Figure [2T] shows 
our Nobeyama NH3 spectra of G35. 03-1-0. 35: while there 
is a narrow NH3(3,3) emission feature that is not well fit 
by the model, the signal-to-noise is insufficient to iden- 
tify it as a candidate maser from the single-dish data. 
This comparison demonstrates that single-dish surveys 
readily miss weak NH3(3,3) masers detected with inter- 
ferometers; sensitive interferometric observations are re- 
quired to assess the prevalence of NH3 masers in MSFRs, 
and their association with other mase r types (see also 
iBroganet aLllMTI : iBrogan et al.l[20Tl . 

4.4. Future Work 

Our analysis of our Nobeyama EGO survey shows that 
the presence of NH3(2,2) and (3,3) emission, H2O masers, 
and Class I and II CH3OH masers are strongly correlated. 
These star formation indicators tend to occur in concert 
(at least on the scales probed by single-dish surveys), 
and identify a (sub)population of EGOs in which central 
MYSO(s) are substantially affecting their environments, 
heating the surrounding gas and exciting maser emis- 
sion. Notably, maser emission and warm dense gas ap- 
pear to pinpoint such sources more effectively than MIR 
indicators such as the "likely" /" possible" classification 
of COB or the presence/ absence of IRDCs. These sources 
are excellent targets for high-resolution foUowup obser- 
vations aimed at understanding the importance of differ- 
ent (proto)stellar feedback mechanisms in massive star 
forming regions, as de monstrated by the SMA, CA RMA, 
and VLA studies o f iCvganowski et al] (|2011al lH) and 
iBrogan et al. I (|2011[ ). These EGOs are also important 
testbeds for proposed maser evolutionary sequences, as 
discussed in more detail below. 

Less clear is the nature of those EGOs detected only 
in NH3(1,1) emission in our survey. An examination 
of their GLIMPSE images suggests they are a hetero- 
geneous group, including both EGOs in IRDCs (e.g. 
G12.02-0.21) and EGOs adjacent to 8 and 24 ^m-bright 
nebulae (e.g. G29.91— 0.81). Some examples of each of 
these MIR source types are detected in H2O maser emis- 
sion, while others are not. The MIR morphologies of 
EGOs without detected H2O masers in our survey are 
similarly heterogeneous, and some H2O maser nonde- 
tections are associated with NH3(2,2) and (3,3) emis- 
sion. Higher-resolution observations are required to lo- 
calize the NH3 and H2O maser emission detected in our 
Nobeyama data with respect to the MIR emission 

We emphasize that high-resolution observations are 
crucial for building an evolutionary sequence for MYSOs, 
and placing EGOs within it. In general, multiple MIR 
sources are present within the Nobeyama beam, and 
detailed studies of EGOs to date reveal mm and cm- 
A multiplicity on '^0.1 pc scales. Furthermore, the 
members of (proto)clusters associated with EGOs ex- 
hibit a range of star formation indicators, suggestive 
of a range of evolutionary states (e.g. .Cvganowski et ahl 
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l2011atlBrogan et al. l[20Tll) . 

EGOs are notably rich in maser emission, and maser 
studies have and continue to provide key insights into 
the nature of EGOs; their copious maser emission hke- 
wise provides opportunities to use EGOs to advance our 
understanding of masers in massive star-forming regions. 
H2O, Class I and II CH3OH, and OH masers are ubiq- 
uitous in regions of massive star formation, and much 
effort has been devoted to placing these different maser 
types into an evolutionary sequence. Of particular inter- 
est is which maser type appears first-and thus pinpoints 
the earliest stages of massive star formation. In most 
proposed sequences. Class I CII3OII masers are identi- 
fied with the earliest stages of MYSO evolution, with the 
youngest sources being those associated only wit h Class 
I CH3OH masers (e.g. EUingsen 2006; EUingse n et al.l 
I2007t [Breen et aLll2010aD . However, recent work suggests 
that Class I CH3OH masers may be excited by shocks 
drive n by expanding HII r egions as well as by outflows 
(e.g. iVoronkov et"aLl [2010l ). such that Class I CH3OH 
may outlast the Class II maser stage and/or arise more 
than on ce during MYSO evolution (e.g.lChen et al. 120111 : 
Voronkov et al.. 201%. IBreen fc EllingsenI (|20TT| ) and 



Casweh fc Breed (|2010D have also recently proposed that 



H2O masers-particularly those with blueshifted high- 
velocity features-may be the earliest signposts of MYSO 
formation, preceding the Class II CH3OH maser stage. 

Statistical comparisons of "Class I only" and "Class II 
only" EGOs based on our data are limited by the small 
sample sizes. It is notable, however, that the NH3(2,2) 
and H2O maser detection rates towards these subsam- 
ples are comparable, particularly considering the small 
number statistics. Likewise, Figures [12] and [18] show 
no clear patterns in their NH3 or H2O maser proper- 
ties that would suggest a trend in evolutionary state. 
The parameter space occupied by Class I-only and Class 
Il-only sources in these plots also largely overlaps with 
that occupied by EGOs associated with both CH3OH 
maser types. Though the comparison is again limited by 
small- number statistics, the difference in the NH3(3,3) 
detection rates (63%/14% for Class I/II-only sources) is 
intriguing, particularly given the association of Class I 
CH3OH and NH3(3,3) masers (e.g. IBrogan eFaLllWl . 

Progress in our understanding of masers as evolution- 
ary indicators for MSF requires identifying candidate 
youngest sources, and studying them in detail. The 
(smaU) samples of EGOs with H20+Class I CH3OH and 
H20-fClass II CH3OH masers identified in our survey 
will be promising targets for such studies, as will the 
samples of H2 0-only, Class I CH30H-only, and Class II 
CH30H-only sources. Sensitive, high-resolution maser 
observations are needed: (1) to localize the maser emis- 
sion, and determine whether or not all maser species are 
associated with the same MYSO and (2) to search for 
weak masers and establish whether maser types unde- 
tected in single-dish surveys are truly absent. The ex- 
panded capabilities of the Karl G. Jansky Very Large 
Array (VLA) are well-suited to such studies. High- 
resolution cm-(sub)mm wavelength line and continuum 
observations will also constrain the properties of com- 
pact cores (temperature, density, mass, chemistry) and 
outflows, allowing maser activity to be correlated with 
other signposts of star formation at the scale of individ- 
ual active sources. 



5. CONCLUSIONS 

We have surveyed all 94 GLIMPSE EGOs visible from 
the northern hemisphere (S > —20°) in H2O maser and 
NH3(1,1), (2,2), and (3,3) emission with the Nobeyama 
45-m telescope. Our results provide strong evidence that 
EGOs, as a population, are associated with dense gas 
and active star formation, and also reveal statistically 
significant variation amongst EGO subsamples: 

• H2O masers, which are associated with outflows 
and require high densities (n(H2)~10^-10^*' cm^'^), 
are detected towards ~68% of EGOs surveyed. 

• The NH3(1,1) detection rate is ^97%, confirming 
that EGOs are associated with dense molecular gas. 

• Two-component models provide a significantly im- 
proved fit for ^23% of our NH3 spectra, indicating 
contributions from both warm inner regions and 
cooler envelopes along the line of sight. 

• H2 O maser emission is strongly correlated with the 
presence of warm, dense g indicated by emis- 
sion in the higher-excitation NH3 transitions. The 
H2O maser detection rate is 81% towards EGOs de- 
tected in NH3(2,2) and/or (3,3) emission, and only 
44% towards EGOs detected only in NH3(1,1). We 
find statistically significant differences in the dis- 
tributions of NH3 column density, kinetic temper- 
ature, and NH3 linewidth for EGOs with/without 
H2O maser detections: EGOs with detected H2O 
masers have greater median N(NH3), Tkm, and 

• H2O maser and NH3(2,2) and (3,3) detection rates 
are higher towards EGOs classified as "likely" out- 
fiow candidates based on their MIR morphology 
than towards EGOs classified as "possible" outfiow 
candidates. However, statistical tests show signif- 
icant differences only in the distributions of the 
NH3(1,1) and (2,2) peaks (Tmb), not in physical 
properties. 

• EGOs associated with IRDCs have higher NH3(2,2) 
and (3,3) detection rates, but a lower H2O maser 
detection rate, than EGOs not associated with 
IRDCs. We find statistically significant differences 
in the distributions of NH3(1,1) peak (Tmb), NH3 
linewidth, and NH3 beam filling fraction for EGOs 
associated/not associated with IRDCs: the median 
NH3(1,1) Tmb is higher, and the median Uv lower, 
for EGOs associated with IRDCs. 

• The H2O maser, NH3(2,2), and NH3(3,3) detec- 
tion rates towards EGOs with both Class I and 
II CH3OH masers are the highest of any EGO 
subsample we consider: 95%, 90% and 81%, re- 
spectively. In contrast, we detect H2O masers 
and the higher-excitation NH3 lines towards only 
33% (H2O), 20% (2,2) and 7% (3,3) of EGOs with 
neither type of CH3OH maser. We find statis- 
tically significant differences in the distributions 
of NH3(1,1) peak temperature (Tmb), NH3 col- 
umn density, and NH3 linewidth for EGOs asso- 
ciated with both types/neither type of CH3OH 
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masers: EGOs associated with both Class I and II 
CH3OII masers have higher median NIl3(l,l) Tmb, 
N(NH3), CTv, and T 

kin ■ 

• While H2O maser detection rates vary across EGO 
subsamples, we find no evidence for statistically 
significant differences in the properties of detected 
H2O masers. 

Our II2O maser and NH3 survey, in combination with 
the 1.1 mm continuum BGPS, provides the necessary 
data to explore connections between H2O maser and 
clump properties: H2O maser spectra, clump-scale Tkin 
measurements from NII3, and clump masses and densi- 
ties from the 1.1 mm dust continuum emission and T^in 
measurements. These combined data show no correlation 
between isotropic II2O maser luminosity and volume- 
averaged clump density. II2 O maser luminosity is weakly 
positively correlated with clump temperature and with 
clump mass, consistent with reported correlations be- 
tween H2O maser luminosity and the bolometric lumi- 
nosity of the driving source. 

We interpret the observed correlations of H2O maser 
and clump properties in terms of the relationship be- 
tween a clump and the massive star(s) it forms. For more 
evolved sources (UC HIIs), studies indicate that higher- 
mass clumps form highe r-mass (and thus more luminous) 
stars (e.g. Johnston et al . 20091. For an actively accret- 
ing MYSO, the available mass reservoir sets the limit on 
its final, stellar mass. The more luminous (and massive) 
an MYSO, the more energy it will impart to its envi- 
rons, and the more it will heat the gas and dust of the 
surrounding clump. 

We find that NH3(2,2) and (3,3) emission, II2O masers, 
and Class I and II CH3OH masers are strongly corre- 
lated, at least on the scales probed by single-dish sur- 
veys. These star formation indicators pinpoint EGOs in 
which the central MYSO(s) are substantially affecting 
their environments, more effectively than MIR indica- 
tors (such as the "likely" /" possible" classification of .COS 
or the presence/absence of IRDCs). We also identify 
small samples of EGOs associated with only one maser 
type; the Il2 0-only and Class I CH30II-only sources are 
candidates for extremely young MYSOs. Constructing 
an evolutionary scheme for MYSOs requires localizing 
maser and dense gas emission at the scale of individ- 
ual (proto)stars. The expanded capabilities of the Karl 
G. Jansky VLA will enable such studies for statistically 
meaningful samples. 
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Table 1 

EGO Sample; Properties from the Literature 



Source Name J2000 Coordinates'^ EGO Catalog^ IRDC?* CH3OH maser?'' 





„ (h m 


S (° ' ") 






Class II 


Class I 


GlO.29-0.13 


18 08 49.3 


-20 05 57 


2 


Y 


Y 


Y 


GlO.34-0.14 


18 09 00.0 


-20 03 35 


2 


Y 


Y 


Y 


Gil. 11-0.11 


18 10 28.3 


-19 22 31 


3 


Y 


Y 


N 


Gil. 92-0.61 


18 13 58.1 


-18 54 17 


1 


Y 


Y 


Y 


G12.02-0.21 


18 12 40.4 


-18 37 11 


1 


Y 


N 


N 


G12.20-0.03 


18 12 23.6 


-18 22 54 


4 


N 


Y 


Y 


G12.42+0.50 


18 10 51.1 


-17 55 50 


4 


N 


N 


Y 


G12.68-0.18 


18 13 54.7 


-18 01 47 


4 


N 


Y 


Y 


G12.91-0.03 


18 13 48.2 


-17 45 39 


1 


Y 


Y 


Y 


G12.91-0.26 


18 14 39.5 


-17 52 00 


5 


N 


Y 


Y 


G14.33-0.64 


18 18 54.4 


-16 47 46 


1 


Y 






G14.63-0.58 


18 19 15.4 


-16 30 07 


1 


Y 


Y 


Y 


G16.58-0.08 


18 21 15.0 


-14 33 02 


3 


Y 


N 


N 


G16.59-0.05 


18 21 09.1 


-14 31 48 


2 


Y 






G16.61-0.24 


18 21 52.7 


-14 35 51 


1 


Y 






G17.96+0.08 


18 23 21.0 


-13 15 11 


4 


N 


N 


N 


G18.67+0.03 


18 24 53.7 


-12 39 20 


1 


N 


Y 


Y 


G18.89-0.47 


18 27 07.9 


-12 41 36 


1 


Y 


Y 


Y 


G19.01-0.03 


18 25 44.8 


-12 22 46 


1 


Y 


Y 


Y 


G19.36-0.03 


18 26 25.8 


-12 03 57 


2 


Y 


Y 


Y 


G19.61-0.12 


18 27 13.6 


-11 53 20 


2 


N 


Y 


N 


G19.61-0.14 


18 27 16.8 


-11 53 51 


4 


N 






G19.88-0.53 


18 29 14.7 


-11 50 23 


1 


Y 


Y 


Y 


G20.24+0.07 


18 27 44.6 


-11 14 54 


4 


N 


Y 


N 


G21. 24+0.19 


18 29 10.2 


-10 18 11 


4 


Y 


N 


N 


G22.04+0.22 


18 30 34.7 


-09 34 47 


1 


Y 


Y 


Y 


G23.01-0.41 


18 34 40.2 


-09 00 38 


1 


N 






G23.82+0.38 


18 33 19.5 


-07 55 37 


4 


N 


- 


Y 


G23.96-0.il 


18 35 22.3 


-08 01 28 


1 


N 


Y 


Y 


G24.00-0.10 


18 35 23.5 


-07 59 32 


1 


Y 


Y 


Y 


G24.11-0.17 


18 35 52.6 


-07 55 17 


4 


Y 


N 


Y 


G24.11-0.18 


18 35 53.0 


-07 55 23 


4 


Y 






G24.17-0.02 


18 35 25.0 


-07 48 15 


1 


Y 


- 


N 


G24.33+0.14 


18 35 08.1 


-07 35 04 


4 


Y 


Y 


Y 


G24.63+0.15 


18 35 40.1 


-07 18 35 


3 


Y 


N 


Y 


G24.94+0.07 


18 36 31.5 


-07 04 16 


1 


N 


Y 


Y 


G25.27-0.43 


18 38 57.0 


-07 00 48 


1 


Y 


Y 


Y 


G25.38-0.15 


18 38 08.1 


-06 46 53 


2 


Y 


- 


Y 


G27.97-0.47 


18 44 03.6 


-04 38 02 


1 


Y 


N 


Y 


G28.28-0.36 


18 44 13.2 


-04 18 04 


2 


N 


Y 


N 


G28.83-0.25 


18 44 51.3 


-03 45 48 


1 


Y 






G28.85-0.23 


18 44 47.5 


-03 44 15 


4 


N 


Y 


N 


G29.84-0.47 


18 47 28.8 


-02 58 03 


3 


Y 






G29.89-0.77 


18 48 37.7 


-03 03 44 


4 


Y 


N 


N 


G29.91-0.81 


18 48 47.6 


-03 03 31 


4 


N 


N 


N 


G29.96-0.79 


18 48 50.0 


-03 00 21 


3 


Y 


N 


N 


G34.26+0.15 


18 53 16.4 


+01 15 07 


5 


N 


- 


Y 


G34.28+0.18 


18 53 15.0 


+01 17 11 


3 


Y 


- 


Y 


G34.39+0.22 


18 53 19.0 


+01 24 08 


2 


Y 


- 


Y 


G34.41+0.24 


18 53 17.9 


+01 25 25 


1 


Y 


- 


Y 


G35.03+0.35 


18 54 00.5 


+02 01 18 


1 


N 






G35.04-0.47 


18 56 58.1 


+01 39 37 


1 


Y 


N 


Y 


G35.13-0.74 


18 58 06.4 


+01 37 01 


1 


N 


- 


Y 


G35.15+0.80 


18 52 36.6 


+02 20 26 


1 


N 


N 


N 


G35.20-0.74 


18 58 12.9 


+01 40 33 


1 


N 


Y 


Y 


G35.68-0.18 


18 57 05.0 


+02 22 00 


1 


Y 


N 


N 


G35.79-0.17 


18 57 16.7 


+02 27 56 


1 


Y 


- 


Y 


G35.83-0.20 


18 57 26.9 


+02 29 00 


4 


Y 






G36.01-0.20 


18 57 45.9 


+02 39 05 


1 


Y 


- 


N 


G37.48-0.10 


19 00 07.0 


+03 59 53 


1 


N 


Y 


N 


G37.55+0.20 


18 59 07.5 


+04 12 31 


5 


N 


- 


N 


G39.10+0.49 


19 00 58.1 


+05 42 44 


1 


N 


Y 


Y 


G39.39— 0.14 


19 03 45.3 


+05 40 43 


4 


N 




Y 


G40.28-0.22 


19 05 41.3 


+06 26 13 


3 


Y 




Y 


G40.28-0.27 


19 05 51.5 


+06 24 39 


1 


Y 


N 


N 


G40.60-0.72 


19 08 03.3 


+06 29 15 


4 


N 


N 


Y 


G43.04-0.45(a) 


19 11 38.9 


+08 46 39 


4 


N 






G43.04-0.45(b) 


19 11 39.1 


+08 46 32 


4 


N 






G44.01-0.03 


19 11 57.2 


+09 50 05 


1 


N 


N 


N 


G45.47+0.05 


19 14 25.6 


+11 09 28 


1 


Y 




N 


G45.47+0.13 


19 14 07.3 


+11 12 16 


4 


N 






G45.50+0.12 


19 14 13.0 


+11 13 30 


4 


N 




N 


G45.80-0.36 


19 16 31.1 


+11 16 11 


3 


N 




Y 


G48.66-0.30 


19 21 48.0 


+13 49 21 


2 


Y 


N 


N 
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Table 1 — Continued 



Source Name J2000 Coordinates'' EGO Catalog'' IRDC?'' CH3OH maser?'' 





„ (h m s) 


5{° 


/ 


") 






Class II 


Class 


G49.07-0.33 


19 22 41.9 


+ 14 


10 


12 


3 


Y 


N 


Y 


G49.27— 0.32 


19 23 02.2 


+ 14 


20 


52 


3 


N 






G49.27— 0.34 


19 23 06.7 


+14 


20 


13 


1 


Y 






G49.42+0.33 


19 20 59.1 


+ 14 


46 


53 


2 


N 


Y 


N 


G49.91+0.37 


19 21 47.5 


+ 15 


14 


26 


4 


N 






G50.36— 0.42 


19 25 32.8 


+15 


15 


38 


3 


Y 






G53.92— 0.07 


19 31 23.0 


+ 18 


33 


00 


3 


N 


N 


N 


G54.il— 0.04 


19 31 40.0 


+ 18 


43 


53 


4 


N 


N 


N 


G54.il— 0.05 


19 31 42.2 


+18 


43 


45 


4 


N 






G54.ll— 0.08 


19 31 48.8 


+ 18 


42 


57 


3 


N 


N 


N 


G54.45+1.01 


19 28 26.4 


+19 


32 


15 


3 


N 


N 


Y 


G54.45+1.02 


19 28 25.7 


+19 


32 


20 


1 


N 






G56. 13+0.22 


19 34 51.5 


+20 


37 


28 


1 


N 






G57.61+0.02 


19 38 40.8 


+21 


49 


35 


4 


N 




N 


G58. 09-0.34 


19 41 03.9 


+22 


03 


39 


1 


Y 






G58. 78+0.64 


19 38 49.6 


+23 


08 


40 


4 


N 






G58. 78+0.65 


19 38 49.2 


+23 


08 


50 


4 


N 






G58. 79+0.63 


19 38 55.3 


+23 


09 


04 


3 


N 






G59. 79+0.63 


19 41 03.1 


+24 


01 


15 


1 


Y 




Y 


G62. 70-0.51 


19 51 51.1 


+25 


57 


40 


3 


N 







^ From lUoa . The table number from[COl is given in the "EGO Catalog" column. Tables 1 and 2 of C08 listed "likely" outflow candidates; 
tables 3 and 4 "possible" outflow candidates. Table 5 sources are those for which only "outflow-only" photometry was presented; we do 
not include them in our analysis of "likely" and "possible" subsamples. 

^ From [UE 111 (only). EGOs in our sample that are not included in lCElll are indicated by ■ ■ • . - indicates a source with a single-dish 
6.7 GHz CH3OH maser detection but no positional information, considered as having "no information" at 6.7 GHz bv lCElJ (see i|3.1.1ll . 



Table 2 

Detection Statistics'' 







NH3 (1,1) 


NH3 (2,2) 


NH3 (3,3)'' 


Water Masers 


Category'^ 




'^detect 


Rate 


^detect 


Rate 


^detect 


Rate 


^detect 


Rate 


Overall 


94 


91 


0.97(0.02) 


59 


0.63(0.05) 


43 


0.46(0.05) 


64 


0.68(0.05) 


IRDC Assoc. 


47 


47 


1.00 


35 


0.74(0.06) 


29 


0.62(0.07) 


29 


0.62(0.07) 


No IRDC Assoc. 


47 


44 


0.94(0.04) 


24 


0.51(0.07) 


14 


0.30(0.07) 


35 


0.74(0.06) 


'Likely' 


48 


47 


0.98(0.02) 


34 


0.71(0.07) 


29 


0.60(0.07) 


36 


0.75(0.06) 


'Possible' 


43 


41 


0.95(0.03) 


22 


0.51(0.08) 


12 


0.28(0.07) 


25 


0.58(0.08) 


'Outflow-only' 


3 


3 


1.00 


3 


1.00 


2 


0.67(0.27) 


3 


1.00 


Detected in: 




















NH3(1,1) only 


32 














14 


0.44(0.09) 


NH3(1,1) and (2,2) 


16 














13 


0.81(0.10) 


NH3(1,1), (2,2) and (3,3) 


43 














35 


0.81(0.06) 


Methanol Maser Associations'^ 


Class I 


41 


41 


1.00 


35 


0.85(0.06) 


31 


0.76(0.07) 


37 


0.90(0.05) 


Class I ND 


28 


25 


0.89(0.06) 


10 


0.36(0.09) 


3 


0.11(0.06) 


13 


0.46(0.09) 


Class II 


28 


27 


0.96(0.04) 


23 


0.82(0.07) 


18 


0.64(0.09) 


24 


0.86(0.07) 


Class II ND 


23 


22 


0.96(0.04) 


9 


0.39(0.10) 


6 


0.26(0.09) 


10 


0.43(0.10) 


Class I Only 


8 


8 


1.00 


6 


0.75(0.15) 


5 


0.63(0.17) 


5 


0.63(0.17) 


Class II Only 


7 


6 


0.86(0.13) 


4 


0.57(0.19) 


1 


0.14(0.13) 


4 


0.57(0.19) 


Class I and II 


21 


21 


1.00 


19 


0.90(0.06) 


17 


0.81(0.09) 


20 


0.95(0.05) 


Neither 


15 


14 


0.93(0.06) 


3 


0.20(0.10) 


1 


0.07(0.06) 


5 


0.33(0.12) 



^ Uncertainties in detection rates calculated using binomial statistics. 
^ Includes only sources detected in NH3(1,1) and (2,2) as well as (3,3), see i|3.1.2l 
IRDC associations and 'likely/possible' designations from IIC09I1 , see also Jl] 

For consistency, all data on Class I and II CH3OH maser associations are taken from Table 1 of ICEIII . see i|3.1.1l 'Class I': all EGOs 
with a Class I maser detection in lCElll . regardless of Class II association (or lack of Class II information). 'Class I ND'; all EGOs listed 
as Class I nondetections in lCElll . regardless of Class II association (or lack of Class II information). 'Class 11': all EGOs listed as Class II 
maser detections in lCEll] . regardless of Class I association. 'Class II ND': all EGOs listed as Class II nondetections in lCElll regardless 
of Class I association. 'Class I Only': EGOs listed as Class I detections and Class II nondetections in lCElll . 'Class II only': EGOs listed 
as Class I nondetections and Class II detections in lCElll . 'Class I and 11': EGOs listed as both Class I and Class II detections in lCElll . 
'Neither': EGOs listed as both Class I nondetections and Class II nondetections in lCElll . 



Table 3 

NH3 Properties: Single Component Fits* 



Source Name 


vlsr 








Distance*^ 


T 


MB 






N(NH3) 






Tex 


Tfci„'' 


Tmb(2,2) = 


Tmb(3,3)'^ 


2 comp?® 


H2O 




(km s-^) 


(km 


s 




(kpc) 




(K) 






(cm-2) 
XIO" 








(K) 


(K) 


(K) 


(K) 




maser? 


GlO.29-0.13 


13.97(0.01) 


1 


18 





01) 


1.58(-l-0. 86,-1. 13) 


1 


94 


0.04) 


4 


76 


22.1(0.6) 





096 


0.002) 


4.52(0.05) 


21.19(0.17) 


1.19(0.04) 


0.61(0.05) 


Y 


Y 


GlO.34-0.14 


12.02(0.02) 


1 


07 





02) 


1.29(-|-0. 92,-1. 23) 


1 


56 


0.06) 


4 


28 


18.2(1.0) 





047 


0.002) 


3.95(0.07) 


28.23(0.38) 


0.94(0.05) 


0.76(0.07) 


Y 


Y 


Gil. 11-0. 11 


29.79(0.02) 





68 





01) 


12.67(-|-0.48,-0.41)f 


1 


79 


0.07) 


5 


15 


11.9(0.6) 





150 


0.006) 


4.45(0.09) 


14.15(0.27) 


0.65(0.06) 


0.32(0.07) 


Y 


N 


Gil. 92-0. 61 


36.11(0.01) 


1 


09 





01) 


3.48(-|-0.44,-0.52) 


3 


02 


0.07) 


6 


40 


31.7(0.8) 





087 


0.001) 


4.79(0.04) 


26.27(0.19) 


1.83(0.06) 


1.14(0.07) 


Y 


Y 


G12.02-0.21 


-3.15(0.06) 


1 


04 





05) 


5. 30(-l-0. 20,-0. 20)8 





46 


0.06) 


5 


18 


13.2(2.1) 





042 


0.004) 


3.16(0.15) 


<12.85 


<0.19 


<0.22 


N 


N 


G12.20-0.03 


51.16(0.06) 


1 


24 





05) 


11.70(+0.31,-0.27)'' 





55 


0.06) 


3 


76 


13.5(2.3) 





032 


0.004) 


3.27(0.13) 


19.56(0.76) 


0.28(0.05) 


<0.25 


N 


Y 


G12.42+0.50 


18.13(0.03) 


1 


16 





04) 


1.87(-|-0.69,-0.85) 





81 


0.06) 


1 


19 


6.0(1.6) 





058 


0.016) 


4.30(0.42) 


29.31(0.60) 


0.55(0.05) 


0.37(0.06) 


N 


Y 


G12.68-0.18 


55.69(0.01) 


1 


18 





01) 


4.46(+0.30,-0.34) 


1 


45 


0.05) 


7 


48 


32.4(0.9) 





051 


0.001) 


3.87(0.04) 


24.95(0.21) 


0.99(0.05) 


0.66(0.05) 


Y 


Y 


G12.91-0.03 


56.79(0.02) 


1 


01 





02) 


4.48(+0.30,-0.34) 


1 


30 


0.05) 


3 


30 


12.6(0.7) 





066 


0.003) 


4.11(0.08) 


23.56(0.31) 


0.72(0.05) 


0.42(0.05) 


Y 


Y 


G12.91-0.26 


37.11(0.02) 


1 


66 





02) 


3.41(+0.43,-0.50) 


2 


00 


0.06) 


4 


66 


31.6(0.9) 





068 


0.002) 


4.30(0.05) 


25.64(0.21) 


1.40(0.06) 


0.79(0.06) 


Y 


Y 


G14.33-0.64 


22.50(0.01) 


1 


23 





01) 


1.12(+0. 13,-0. 13)'' 


2 


35 


0.06) 


3 


58 


20.9(0.5) 





105 


0.002) 


5.11(0.06) 


25.26(0.17) 


1.49(0.06) 


0.79(0.06) 


Y 


Y 


G14.63-0.58 


18.64(0.01) 





87 





01) 


1.73(+0. 63,-0. 76) 


1 


65 


0.06) 


3 


29 


11.5(0.6) 





109 


0.005) 


4.71(0.10) 


20.76(0.32) 


0.95(0.07) 


0.45(0.07) 


N 


Y 


G16.58-0.08 


41.06(0.02) 





70 





02) 


3.24(+0.39,-0.44) 





65 


0.04) 


5 


28 


9.8(0.8) 





060 


0.004) 


3.38(0.12) 


13.56(0.47) 


0.25(0.04) 


<0.17 


N 


N 


G16.59-0.05 


59.84(0.02) 


1 


22 





02) 


4.20(-l-0.30,-0.33) 





91 


0.04) 


2 


10 


9.2(0.8) 





078 


0.006) 


4.13(0.12) 


20.51(0.38) 


0.46(0.04) 


0.33(0.05) 


N 


Y 


G16.61-0.24 


43.61(0.03) 





88 





03) 


3.39(-l-0.37,-0.42) 





69 


0.04) 


3 


23 


8.4(0.8) 





049 


0.004) 


3.48(0.11) 


18.07(0.55) 


0.34(0.04) 


0.19(0.04) 


N 


N 


G17.96+0.08 


22.75(0.05) 





70 





05) 


1.87(4-0.53,-0.61) 





31 


0.04) 





04 


0.4(0.0) 


1 


000 


0.000) 


20.13(0.00) 


<20.13 


<0.17 


<0.17 


N 


Y 


G18.67+0.03 


80.09(0.04) 


1 


23 





03) 


10.81(+0.24,-0.23)' 





58 


0.05) 


5 


34 


20.1(1.7) 





027 


0.002) 


3.26(0.08) 


22.38(0.56) 


0.36(0.05) 


0.21(0.05) 


N 


Y 


G18.89-0.47 


66.17(0.01) 


1 


32 





01) 


4.28(4-0.28,-0.31) 


2 


25 


0.05) 


4 


97 


30.1(0.7) 





076 


0.001) 


4.68(0.04) 


28.24(0.19) 


1.50(0.06) 


1.23(0.07) 


Y 


Y 


G19.01-0.03 


59.92(0.02) 





86 





02) 


11. 54(4-0. 29,-0. 27)'''' 


1 


09 


0.06) 


4 


27 


12.8(1.0) 





047 


0.003) 


3.72(0.08) 


23.79(0.43) 


0.67(0.05) 


0.37(0.05) 


Y 


Y 


G19.36-0.03 


26.79(0.02) 





98 





02) 


2.09(4-0.48,-0.54) 


1 


64 


0.06) 


4 


15 


16.2(0.7) 





068 


0.002) 


4.24(0.07) 


24.90(0.31) 


0.93(0.07) 


0.62(0.07) 


Y 


N 


G19.61-0.12 


57.29(0.08) 


1 


30 





08) 


3.84(4-0.31,-0.34) 





31 


0.05) 


2 


74 


10.9(3.3) 





015 


0.003) 


3.07(0.15) 


25.63(1.21) 


0.24(0.05) 


<0.21 


N 


Y 


G19.61-0.14 


58.38(0.13) 


2 


14 





11) 


3.89(4-0.31,-0.34) 





35 


0.05) 


4 


42 


28.9(5.0) 





Oil 


0.001) 


2.99(0.10) 


26.53(0.92) 


0.22(0.05) 


<0.20 


N 


Y 


G19.88-0.53 


43.50(0.04) 


1 


10 





03) 


3.12(4-0.37,-0.41) 





94 


0.07) 


3 


86 


15.1(1.5) 





039 


0.003) 


3.63(0.10) 


25.72(0.59) 


0.55(0.08) 


0.37(0.08) 


N 


Y 


G20.24-I-0.07 


70.71(0.04) 





74 





03) 


4.38(4-0.28,-0.30) 





78 


0.08) 


/t 
4 


yo 


11.2(1.5) 





047 


0.005) 


3.52(0.17) 


19.23(0.87) 


0.51(0.08) 


<0.33 


N 


Y 


G21. 24-1-0. 19 


25.71(0.03) 





65 





02) 


1.92(4-0.47,-0.53) 





61 


0.05) 


/t 

4 


00 


7.8(0.8) 





058 


0.005) 


3.38(0.14) 


<13.76 


<0.19 


<0.19 


N 


N 


G22. 04-1-0. 22 


51.34(0.02) 





80 


g 


02) 


3.40(4-0.34,-0.37) 





73 


0.04) 




18 


11.5(1.0) 





029 


0.002) 


3.44(0.07) 


26.71(0.49) 


0.49(0.05) 


0.37(0.04) 


Y 


Y 


G23.01-0.41 


77.18(0.02) 


1 


49 





02) 


4.59(+0.38,-0.33)'' 


1 


98 


0.06) 


cr 



( 


34.0(1.0) 





070 


0.001) 


4.25(0.04) 


24.37(0.21) 


1.19(0.06) 


0.66(0.06) 


Y 


Y 


G23.82-l-0.38j 


75.99(0.14) 


1 


48 





13) 


4.43(+0.29,-0.30) 





15 


0.03) 


4 


10 


22.4(6.3) 


0.002(<0.001) 


2.80(0.36) 


<44.20 


<0.10 


0.11(0.03) 


N 


Y 


G23.96-0.il 


73.14(0.02) 


1 


44 





02) 


10.76(+0.28,-0.27)'' 


1 


02 


0.04) 


5 


86 


29.7(1.2) 





033 


0.001) 


3.54(0.04) 


26.78(0.29) 


0.58(0.05) 


0.50(0.05) 


Y 


Y 


G24.00-0.10 


71.70(0.04) 


1 


29 





03) 


4.25(4-0.29,-0.31) 





62 


0.04) 


3 


02 


11.9(1.1) 





043 


0.003) 


3.47(0.10) 


19.60(0.53) 


0.35(0.05) 


<0.19 


N 


Y 


G24.11-0.17 


80.81(0.02) 





81 





02) 


4.62(4-0.28,-0.29) 


1 


22 


0.05) 


3 


05 


8.3(0.5) 





105 


0.005) 


4.18(0.10) 


16.49(0.33) 


0.52(0.05) 


0.26(0.05) 


Y 


N 


G24.11-0.18 


80.79(0.02) 





85 





02) 


4.62(4-0.28,-0.29) 


1 


18 


0.04) 


2 


73 


7.9(0.5) 





106 


0.006) 


4.23(0.10) 


16.77(0.33) 


0.55(0.05) 


0.19(0.05) 


N 


N 


G24.17-0.02 


100.32(0.05) 





56 





05) 


5.41(4-0.29,-0.28) 





17 


0.02) 





01 


0.2(0.0) 


1 


000 


0.000) 


25.40(0.00) 


25.40(1.62) 


0.11(0.03) 


<0.10 


N 


Y 


G24.33-I-0.14 


113.69(0.01) 


1 


03 





01) 


9.16(4-0.29,-0.32)'' 


1 


61 


0.06) 


7 


10 


30.0(0.8) 





050 


0.001) 


4.04(0.04) 


28.82(0.26) 


1.08(0.06) 


1.02(0.06) 


Y 


Y 


G24.63-I-0.15 


53.25(0.02) 





78 





02) 


3.37(4-0.34,-0.36) 


1 


08 


0.05) 


4 


08 


10.1(0.6) 





098 


0.005) 


3.98(0.10) 


15.39(0.38) 


0.47(0.06) 


0.34(0.05) 


Y 


Y 


G24.94-I-0.07 


42.04(0.08) 


1 


28 





07) 


2.77(4-0.37,-0.40) 





33 


0.05) 


3 


88 


16.1(3.4) 





Oil 


0.001) 


3.02(0.12) 


29.73(1.16) 


0.24(0.06) 


0.24(0.05) 


N 


Y 


G25.27-0.43 


59.40(0.03) 





58 





03) 


11.22(4-0.31,-0.30)'' 





88 


0.07) 


4 


43 


7.3(0.9) 





079 


0.008) 


3.64(0.19) 


<14.15 


<0.30 


<0.30 


N 


Y 


G25.38-0.15 


95.47(0.02) 


1 


00 





02) 


5.21(+0.30,-0.29) 


1 


26 


0.07) 


4 


86 


18.6(1.1) 





047 


0.002) 


3.88(0.08) 


26.89(0.44) 


0.83(0.08) 


0.61(0.08) 


N 


Y 


G27.97-0.47 


45.36(0.03) 





54 





03) 


2.86(+0.37,-0.39) 





56 


0.05) 


4 


54 


6.6(1.0) 





046 


0.005) 


3.32(0.19) 


15.60(0.88) 


0.30(0.06) 


0.29(0.06) 


N 


Y 


G28.28-0.36 


49.14(0.02) 





88 





02) 


3.05(+0.36,-0.38) 


1 


21 


0.07) 


6 


28 


17.6(1.0) 





084 


0.004) 


3.91(0.10) 


16.79(0.42) 


0.79(0.08) 


<0.31 


Y 


N 


G28.83-0.25 


86.97(0.03) 


1 


05 





03) 


4.88(+0.34,-0.33) 


1 


09 


0.07) 


3 


69 


14.8(1.3) 





040 


0.003) 


3.77(0.09) 


28.27(0.50) 


0.84(0.08) 


0.67(0.07) 


N 


N 


G28.85-0.23 


96.62(0.06) 





97 





06) 


5.35(+0.38,-0.35) 





38 


0.05) 


3 


86 


9.6(1.8) 





030 


0.004) 


3.12(0.21) 


<15.59 


<0.21 


<0.22 


N 


N 


G29.84-0.47 


68.90(0.06) 





72 





06) 


4.02(+0.34,-0.34) 





34 


0.05) 


3 


45 


6.3(1.9) 





031 


0.006) 


3.11(0.28) 


<14.88 


<0.21 


<0.21 


N 


N 


G29.89-0.77 


83.90(0.04) 





97 





04) 


4.75(4-0.36,-0.34) 





64 


0.06) 


3 


64 


9.9(1.4) 





054 


0.006) 


3.45(0.17) 


<15.79 


<0.28 


<0.27 


N 


N 


G29.91-0.81 


84.33(0.06) 





72 


u 


06) 


4.77(4-0.36,-0.34) 





56 


0.07) 





4 ( 


4.5(1.8) 





063 


0.018) 


3.31(0.43) 


<11.91 


<0.27 


<0.28 


N 


N 


G29.96-0.79 


85.13(0.02) 





72 





02) 


4.81(4-0.36,-0.35) 





54 


0.03) 


1 


71 


3.7(0.6) 





066 


0.009) 


3.66(0.13) 


16.71(0.47) 


0.19(0.03) 


0.11(0.03) 


N 


N 


G34.26-I-0.15 


58.56(0.01) 


1 


77 





01) 


3.52(4-0.39,-0.39) 


2 


59 


0.06) 


3 


30 


33.9(0.6) 





078 


0.001) 


5.34(0.05) 


36.20(0.20) 


2.21(0.06) 


2.02(0.06) 


Y 


Y 


G34.28-I-0.18 


56.09(0.03) 





89 





03) 


3.38(4-0.38,-0.39) 





70 


0.05) 


2 


48 


6.9(1.1) 





059 


0.007) 


3.69(0.15) 


18.69(0.65) 


0.34(0.05) 


<0.19 


N 


Y 


G34. 39-1-0. 22 


57.65(0.01) 


1 


14 





01) 


1.56(4-0.12,-0.11)'' 


1 


88 


0.05) 


4 


38 


20.3(0.6) 





088 


0.002) 


4.52(0.06) 


22.89(0.22) 


1.12(0.06) 


0.53(0.06) 


Y 


Y 


G34.41+0.24 


57.92(0.01) 


1 


18 





01) 


1.56(+0. 12,-0. 11)'' 


2 


62 


0.05) 


4 


64 


25.7(0.5) 





092 


0.001) 


4.93(0.04) 


26.49(0.16) 


1.47(0.05) 


1.07(0.06) 


Y 


Y 


G35.03-t-0.35 


53.37(0.06) 


1 


19 





06) 


3.24(+0.39,-0.39) 





41 


0.05) 


2 


15 


8.8(2.4) 





019 


0.004) 


3.25(0.13) 


29.54(0.92) 


0.33(0.05) 


0.27(0.05) 


N 


Y 


G35.04-0.47 


51.15(0.05) 





65 





05) 


3.11(+0.39,-0.40) 





44 


0.05) 


1 


81 


3.2(1.4) 





058 


0.018) 


3.42(0.31) 


<14.59 


<0.20 


<0.21 


N 


N 


G35.13-0.74 


33.79(0.02) 


1 


82 





01) 


2.11(+0.41,-0.43) 


2 


39 


0.06) 


3 


35 


28.0(0.7) 





116 


0.003) 


5.14(0.06) 


23.34(0.17) 


1.22(0.06) 


0.77(0.06) 


Y 


Y 


G35.15+0.80 
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Table 3 — Continued 



Source Name 


vlsr 




Distance*^ 


TmbCI,!)" T(i 1) N(NH3) 


Vff 






Tmb(2,2)<^ TmbCS.S)" 2 comp?" H2O 




(km s"^) 


(km s-^) 


(kpc) 


(K) (cm-2) 
xio" 




(K) 


(K) 


(K) (K) mascr? 



^ Uncertainties are given in parentheses. For kinematic distances, the uncertainties are based on the prescription of iReid et al.| l i200g|[) (see also j|3.2. 1| >: for maser parallax distances, the 
uncertainties are taken from the cited reference. For the NH3(1,1), (2,2), and (3,3) peak temperatures, the quoted uncertainty is the Icr rms. For all other quantities, the uncertainties are 
estimated from the model optimization, and uncertainties of 0.00 indicate the case Tfci^ ^Te:c (sec jfm . 
^ Near kinematic distance estimated using the NH3 vlsr, except as otherwise noted. See also ^3.2.11 
^ Peak temperature of the NH3 emission on the Tmb scale. All upper limits are 4(7. 
^ Tfci^ is indicated as an upper limit if NH3(2,2) emission is not detected at >4(T. 

^ Indicates whether a two component model was fit (" ^2.2| >. If 'Y', the two-component model results arc listed in Table [4] 

Associated with a 6.7 GHz CH3OH maser assigned the far distance by lGreen fe McClure-Griffithj 120111) . Except for G12.20-0.03 and G45.80-0.36, all far distance assignments 
are "b" classifications in their scheme (sec also j|3.2. It . We adopt the far kinematic distance estimated from the NH3 vlsr- 

° The longitude and velocity of this source indicate that it is likely in the near 3 kpc arm (sec for example Fig. 1 of I Green" et al .112003 ) . Following iGreen Si McClure-GriERthg |[20Tll ). 
wc place this source on a circle of radius 3.4 kpc around the Galactic Center, and adopt a distance uncertainty of ±0.2 kpc. 

^ Maser parallax distance. References: G14.33-0.64. ISato et al.l ||2010I ). G23.01-0.41. IBrunthaler et al~l 1120091 ). G34.39+0.22 and G34.41+0.24: IKuravama et all 120111) . 
G35.20-0.74. IZhang eFall {200^) . 
It is unclear if this source is at the near or the far distance (e.g. ICyganowski et ani2011al : IGreen & McClure-Griffithj[20Tll) . 

Source that meets 4(T detection criterion for NH3(1,1) and (3,3), but not (2,2); hence, T^^^ is treated as an upper limit as for other (2,2) nondetections. See also ^3. 1.2| 
NH3 nondetection. 049.42+0.33: distance estimated using H'^^C0+ velocity from lC09l . G53.92-0.07: distance estimated using H2O maser peak velocity. G57.61+0.02: distance 
estimated from velocity of weak (3.9cr) NH3(1,1) emission below our 4cr detection threshold. The distance for G57. 61+0.02 is included here and in Figure[4]for completeness, but 
this source is otherwise excluded from our analysis. See also ^3.2.11 

^ This temperature limit is likely excessively low because the criteria of least-squares fitting are failing and the error distributions are non-Gaussian. 
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Table 4 

NH3 Properties: Two Component Fits^ 



Source Name 




^LSR 


CTv 


Tea; 


N(NH3) 




(K) 


(km s-l) 


(km s-l) 


(K) 


(cm~^) 



GlO.29-0.13 
GlO.34-0.14 
Gil. 11-0. 11 
Gll.92-0.61 
G12.68-0.18 
G12.91-0.03 
G12.91-0.26 
G18.89-0.47 
G19.01-0.03 
G22.04+0.22 
G23.01-0.41 
G23.96-0.il 
G24.11-0.17 
G24.33-f0.14 
G24.63-I-0.15 
G28.28-0.36 
G34.26-f0.15 
G34.39-f0.22 
G34.4H-0.24 
G35.13-0.74 
G35.20-0.74 



15.93(0.32) 
32.91(0.98) 
15.62(0.78) 
41.51(2.42) 
13.90(0.58) 
14.07(0.64) 
13.55(0.20) 
54.06(1.38) 
14.19(0.48) 
31.47(0.78) 
16.77(0.75) 
32.50(1.44) 
18.86(0.53) 
28.48(0.40) 
13.46(0.38) 
52.40(1.92) 
14.74(0.50) 
50.19(4.44) 
15.79(0.68) 
67.75(12.78) 
14.25(0.44) 
30.00(0.49) 
12.06(0.89) 
34.03(1.08) 
12.71(0.74) 
29.26(1.22) 
27.25(0.38) 
34.59(1.53) 
12.27(0.59) 
40.49(4.15) 
10.61(1.11) 
24.15(0.93) 
34.72(0.35) 
43.70(1.33) 
15.57(0.40) 
40.24(1.85) 
14.16(0.31) 
61.20(3.27) 
4.22(0.20) 
25.59(0.20) 
15.19(0.89) 
60.37(9.64) 



14.16(0.02) 
13.31(0.06) 
12.05(0.03) 
12.05(0.05) 
29.40(0.04) 
30.33(0.14) 
36.10(0.01) 
35.89(0.06) 
56.42(0.04) 
54.99(0.04) 
56.77(0.02) 
56.90(0.07) 
37.60(0.02) 
36.48(0.05) 
66.34(0.02) 
65.67(0.05) 
59.95(0.02) 
59.26(0.13) 
51.44(0.02) 
50.77(0.16) 
75.81(0.13) 
77.36(0.02) 
73.04(0.04) 
73.20(0.05) 
80.83(0.02) 
80.76(0.07) 
113.72(0.01) 
113.67(0.11) 
53.30(0.02) 
52.93(0.14) 
49.58(0.03) 
48.47(0.12) 
58.54(0.01) 
58.65(0.07) 
57.72(0.02) 
57.37(0.07) 
58.01(0.01) 
57.60(0.05) 
32.35(0.02) 
34.09(0.02) 
34.21(0.02) 
33.98(0.05) 



0.89(0.02) 
1.73(0.04) 
0.70(0.04) 
1.43(0.05) 
0.45(0.03) 
0.65(0.06) 
0.73(0.01) 
3.43(0.07) 
0.90(0.02) 
1.08(0.03) 
0.74(0.03) 
1.56(0.06) 
0.62(0.02) 
2.13(0.04) 
1.03(0.02) 
2.46(0.05) 
0.62(0.02) 
2.47(0.18) 
0.64(0.03) 
2.35(0.25) 
3.06(0.11) 
0.84(0.02) 
0.89(0.04) 
1.81(0.05) 
0.56(0.02) 
1.32(0.06) 
0.78(0.02) 
2.45(0.10) 
0.65(0.02) 
1.75(0.13) 
0.46(0.04) 
1.11(0.06) 
1.34(0.03) 
2.96(0.07) 
0.87(0.02) 
1.99(0.07) 
1.00(0.01) 
2.36(0.05) 
0.50(0.02) 
1.79(0.02) 
1.28(0.03) 
1.97(0.06) 



4.19(0.06) 
3.23(0.12) 
4.04(0.17) 
3.32(0.23) 
3.70(0.25) 
3.72(0.22) 
5.16(0.07) 
54.06(0.00) 
3.55(0.11) 
3.33(0.10) 
4.35(0.20) 
3.02(0.16) 
4.21(0.12) 
3.47(0.06) 
4.44(0.09) 
52.40(0.00) 
3.65(0.12) 
3.08(0.35) 
3.43(0.14) 
3.07(0.74) 
3.45(0.09) 
3.81(0.08) 
3.43(0.22) 
3.15(0.12) 
4.14(0.21) 
2.97(0.14) 
3.90(0.05) 
2.98(0.13) 
3.81(0.17) 
2.90(0.31) 
3.29(0.28) 
3.73(0.22) 
5.05(0.10) 
3.32(0.12) 
4.27(0.09) 
3.15(0.17) 
5.00(0.08) 
61.20(0.00) 
4.22(0.00) 
4.75(0.06) 
4.46(0.14) 
60.37(0.00) 



15.6(0.7) 
30.2(2.8) 
6.5(1.1) 
21.3(2.7) 
10.6(1.1) 
8.8(1.2) 
19.3(0.6) 
3.4(0.1) 
16.5(1.0) 
36.7(2.0) 
5.3(0.8) 
36.4(6.3) 
8.9(0.7) 
51.6(2.8) 
20.3(0.6) 
3.0(0.1) 
10.3(0.9) 
13.8(7.1) 
8.5(0.8) 
7.4(7.9) 
37.1(2.7) 
23.0(1.2) 
9.7(1.4) 
39.4(3.9) 
4.4(0.6) 
33.3(6.5) 
20.8(1.0) 
78.4(9.7) 
8.2(0.7) 
31.4(8.5) 
27.1(7.0) 
9.5(2.0) 
19.7(1.0) 
63.4(5.2) 
13.6(0.7) 
33.0(4.1) 
17.4(0.4) 
2.3(0.1) 
3.5(0.5) 
28.2(0.9) 
13.1(0.5) 
1.4(0.3) 



Uncertainties estimated from the model optimization are given in parentheses: values of 0.00 indicate cases where the model is poorly constrained 
(sec JO}. 



Table 5 

EGO Subsamples with Statistically Significant Differences in NH3 Properties 



Subs ample 


Property 


K-S significj 


Likely/Possible 


NH3(1,1) peak (Tmb) 


8.4E-03 


Likely /Possible 


!L\tis[z,z) peak (ImbJ 


V at? HQ 
( .OHj-Uo 


iJX-UL^/no IrtiJC' 


JNhl3(l,ij peak (ImbJ 


1.0±!j-U0 


IRDC/no IRDC 




8.2E-04 


IRDC/no IRDC 


Vff 


1.7E-03 


H2O maser detections/ nondetections 




O.Oi1j-U<5 


H2O maser detections/nondetections 


IN (^IN rl3 ) 


1 fit? r\A 


H2O maser detections/nondetections 


kin 




Class I/Class I ND 


NH3(1,1) peak (Tmb) 


3.0E-03 


Class I/Class I ND 




5.1E-03 


i^iass i/oiass i inu 


IN [i\ Jrl3 ) 




L^iass i/oiass i inu 


kin 


K 1 TP HQ 

O.lllj-Uo 


Class Il/Class II ND 


NH3(1,1) peak (Tmb) 


3.2E-03 


Class Il/Class II ND 


CTv 


4.0E-03 


Class Il/Class II ND 


N(NH3) 


2.7E-06 


Class Il/Class II ND 




6.8E-03 


Class Il/Class II ND 


NH3(2,2) peak (Tmb) 


3.9E-03 


Class Il/Class II ND 


NH3(3,3) peak (Tmb) 


3.5E-03 


Both Class I & II/Neither 


NH3(1,1) peak (Tmb) 


8.9E-04 


Both Class I & II/Neither 




1.6E-03 


Both Class I & II/Neither 


N(NH3) 


3.9E-05 
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Table 6 

H2O Maser Properties: Detections 



Source Name 






Vmax 


^pcak 


VrangG 




Speak 




Liso(H2 0) 


Notes^ 




Jy 


km S~^ 


km s~^ 


km 


km s~^ 


km 


Jy 


Jy km s"-"^ 


L© 




1 n on n 1 Q 


0.08 


5.3 


16.1 


12.3 


10.8 


14.0 


0.7 


1.548 


QnTP no 

o.oyii/-Uo 


I3T71 1 

rSlL/11 


CjiU.v>4-U. 14 


0.11 


-13.1 


49.8 


19.8 


62.8 


12.0 


11.1 


a A C\A A 

D4.U44 


Z.40i1;-UD 


"Dm 1 
rilL/ll 


1 1 no n 1 


0. 14 


18.2 


43.8 


39.8 


25.6 


36.1 


53.1 


-\ AC Q /I /I 

140.044 


A CidTT' nc; 
4.UDii(-Uo 


rlGyb,l3r!jll 


1 on 00 

Glz.2U-U.Uo 


0.11 


47.3 


47.6 


47.3 


0.3 


51.2 


0.5 


0.242 


7.bzlir-U7 


BEll 


/^io /10 1 n c;n 
L:j1z.4z+U.oU 


0.11 


4.5 


11.2 


5.8 


6.7 


18.1 


2.9 


7.828 


n QnTT' n'7 

b.oUrLr-U/ 




Cj1z.Do-U.1o 


0.11 


8.3 


109.7 


59.5 


101.4 


55.7 


629.0 


1 Q0 1 Qon 
1081.0ZO 


n QOTT* C\A 

b.ozii(-U4 


T3T71 1 


1 ni n HQ 
Cjlz.yl-U.Uo 


0.11 


13.2 


16.9 


16.7 


3.8 


56.8 


0.7 


1.327 


a. 1 Q TP nv 




G12.91-0.26 


0.14 


40.4 


49.3 


40.7 


8.9 


37.1 


1.0 


1.010 


2.70E-07 


* 


G14.33-0.64 


0.13 


14.3 


32.3 


27.2 


18.1 


22.5 


35.8 


108.120 


3.12E-06 


SIO 


G14.63-U.58 


0.14 


21.8 


22.6 


22.3 


0.8 


18.6 


1.7 


1.452 


l.OOE-07 




G16. 59-0.05 


0.10 


50.1 


72.5 


57.6 


22.4 


59.8 


37.8 


167.848 


6.81E-05 


BEll 


G17.96+0.08 


0.10 


21.0 


31.8 


28.8 


10.8 


22.7 


29.5 


43.348 


3.49E-06 




G18.67+0.03 


0.10 


89.3 


92.3 


90.7 


3.0 


80.1 


5.7 


8.142 


2.19E-05 




G18. 89-0.47 


0.13 


61.3 


69.9 


62.6 


8.6 


66.2 


4.7 


13.412 


5.65E-06 




G19.01-0.03 


0.10 


49.5 


56.8 


54.9 


7.3 


59.9 


1.0 


2.347 


7.19E-06 




G19.61-0.12 


0.12 


57.2 


58.6 


57.8 


1.3 


57.3 


2.3 


2.470 


8.38E-07 




G19.61-0.14 


0.12 


57.2 


58.3 


57.8 


1.1 


58.4 


1.5 


1.472 


5.12E-07 




G19. 88-0.53 


0.18 


20.6 


46.5 


43.8 


25.9 


43.5 


28.6 


119.711 


2.68E-05 


B02 


G20. 24+0.07 


0.19 


69.9 


71.8 


71.3 


1.9 


70.7 


1.5 


2.560 


1.13E-06 




G22.04+0.22 


0.10 


43.2 


54.0 


53.5 


10.8 


51.3 


5.8 


8.327 


2.21E-06 




G23.01-0.41 


0.14 


61.0 


82.3 


69.7 


21.3 


77.2 


52.5 


268.814 


1.30E-04 


FC99 


Gzo.oz+U.oo 


0.06 


72.5 


76.3 


74.4 


3.8 


76.0 


1.5 


2.927 


1 QOTP C\Cl 

Lozt-Ub 




OQ aa nil 
Gzo.yb-U.ll 


0.11 


-165.2 


73.4 


72.6 


238.6 


73.1 


0.9 


3.629 


y.bbii(-Ub 


T3 1 1 

rSll 


A nn n 1 n 
Gz4.UU-U.lU 


0.11 


61.8 


64.8 


64.0 


3.0 


71.7 


1.4 


2.020 


Q QOTT* nT 




A 1 '7 n no 
Gz4. 1 /-U.Uz 


0.06 


92.5 


94.1 


93.3 


1.6 


100.3 


1.6 


1.758 


1 1 oTT* nc 
1. lorlj-Ub 




A QQ 1 n 1/1 

Gz4.oo-hU.14 


0.17 


38.9 


131.1 


110.9 


92.2 


113.7 


12.5 


47.968 


n o/^TT* nc; 


GGll 


A 1 n ic; 


0. 14 


28.6 


31.9 


28.9 


3.2 


53.3 


0.9 


1.745 


A c;/^T7' nT 
4.0bll;-U/ 




G24.y4-hU.U7 


0.13 


30.0 


49.9 


42.9 


20.0 


42.0 


96.5 


125.347 


1 TT* n c 




o"? n /I Q 


0.17 


53.8 


56.5 


54.6 


2.7 


59.4 


2.3 


5.169 


1 ccniT' nc; 
1.0Uii(-U0 




/^oc; QQ n ic; 
Gzo.oo-U.lO 


0.18 


54.6 


97.2 


96.9 


42.6 


95.5 


1.4 


3.526 


onTP r\ci 
z.zUll/-Ub 




c^^i n'7 n /I 


0.13 


46.8 


74.8 


72.1 


28.0 


45.4 


13.7 


61.551 


1. IbUz-Uo 


WUb 


G34.26+0.15 


0.13 


-46.0 


120.1 


60.8 


166.1 


58.6 


182.8 


1167.865 


3.33E-04 




CTo4.za+U.io 


0.10 


60.5 


63.7 


60.8 


3.2 


56.1 


0.8 


0.776 


n /I TP nT 

Z.U4ll;-U / 




Q /I Qn 1 n 00 


0. 14 


30.8 


65.1 


61.6 


34.2 


57.7 


20.8 


97.743 


c: A TIP r\Ci 

o.4^1i/-Ub 


WUb 


Cjv)4.4i+U.z4 


0.13 


8.2 


103.1 


56.5 


94.9 


57.9 


219.2 


OOQ ^ 0'7 

823.427 


4.61E-05 


W06 


Ctoo.Uo+U.oo 


0.11 


36.8 


60.0 


44.6 


23.2 


53.4 


31.1 


62.190 


1 c^nTP nc; 
I.0UII/-U0 


r Gyy 


Ooo.io-U. /4 


0. 14 


-20.9 


85.3 


36.0 


106.2 


33.8 


811.8 


ncif\t^ 1 n /I 
2690.104 


TtlTP C\A 

z. / bll/-U4 




Ooo. iO+U.oU 


n 1 


DD.Z 


vn 1 

/y.i 




ion 
iz.y 


7/1 Q 


1 Q 
Iz.O 


1 Q ono 
io.zUz 


0.0 1 Hz-Ub 




Q on n v/i 


n 1 n 


01 

zi.y 


QQ n 

00. y 


Q1 1 

ol.i 


170 


Q/1 

o4.z 


1 /I 

lz.4 


Qfi Qn/1 


4.UUllj-Ub 


r Gyy 


tjoo.Do-U.lo 


nil 
U.ii 


ov n 
z ( .y 


Qn n 
oU.y 


on Q 

zy.o 


Q n 
o.U 


OQ 1 

zo. i 


0. / 


Q /I 


fi nQT7 nT 
b.Uoiij-U / 




Ooo. / y-u. i / 


n 1 c 


/ .0 


D4.D 


ftQ Q 
00.0 


T n 


^;i Q 


D.4 


1 n ni 
ly .Ulz 


ncT ^^; 
b.UoHz-Ub 


"Dl 1 

rSll 


/^Qv /IC n in 
00/ .40-U.1U 


n n'7 
U.U r 


/inn 
49. U 


Di. / 


f^Q n 
Do.y 


1 v 
iz. 1 


Do.O 


A Q 
4.0 


n n/if? 

y.y4o 


nc;i7 r\ci 
z.yoll/-Ub 




G37.55+0.20 


0.17 


83.5 


89.7 


84.8 


6.2 


85.4 


2.9 


9.951 


7.36E-06 


* 


G39. 10+0.49 


0.08 


23.1 


24.7 


23.9 


1.6 


23.4 


1.8 


2.180 


1.08E-07 


Bll 


G39.39-0.14 


0.08 


62.5 


65.7 


64.9 


3.2 


66.3 


0.4 


0.613 


2.51E-07 


B02 


G40.28-0.22 


0.07 


63.9 


97.3 


72.8 


33.4 


73.6 


59.3 


131.506 


7.50E-05 




G43. 04-0. 45(a) 


0.09 


45.7 


61.9 


48.4 


16.2 


57.7 


2.0 


8.053 


2.83E-06 




G43. 04-0. 45(b) 


0.14 


-11.0 


58.0 


57.2 


69.0 


57.4 


1.3 


6.395 


2.23E-06 




G45.47+0.05 


0.09 


47.9 


74.0 


50.8 


26.2 


61.5 


4.0 


16.435 


8.60E-06 


FC99 


G45.47+0.13 


0.12 


56.6 


67.7 


57.1 


11.1 


60.9 


1.6 


1.946 


9.85E-07 


FC99 


G45.50+0.12 


0.15 


57.7 


74.1 


58.8 


16.4 


61.0 


3.7 


7.566 


3.84E-06 




G45.80-0.36 


0.12 


55.3 


64.4 


63.1 


9.2 


58.8 


8.4 


11.527 


1.28E-05 




G49. 07-0.33 


0.14 


57.7 


58.5 


58.0 


0.8 


60.8 


0.9 


0.837 


5.82E-07 




G49.27-0.34 


0.18 


69.3 


70.1 


69.6 


0.8 


67.9 


1.8 


1.416 


9.78E-07 




G49. 42+0.33 


0.08 


-182.6 


-25.4 


-28.7 


157.2 




2.9 


7.018 


2.33E-05 




G49.91+0.37 


0.11 


5.5 


5.8 


5.5 


0.3 


8.4 


0.6 


0.317 


2.05E-09 




G50.36-0.42 


0.05 


45.2 


45.5 


45.2 


0.3 


39.1 


0.2 


0.128 


2.69E-08 




G53.92-0.07 


0.05 


45.8 


47.1 


46.9 


1.3 




0.3 


0.430 


2.42E-07 




G54.11-0.08 


0.10 


28.2 


48.2 


36.0 


20.0 


39.7 


12.1 


70.951 


2.42E-05 




G54.45+1.01 


0.16 


-33.1 


37.6 


-5.6 


70.6 


35.6 


3.7 


14.909 


3.56E-06 




G54.45+1.02 


0.17 


-33.1 


71.0 


-5.6 


104.1 


35.8 


5.3 


25.161 


6.08E-06 




G58.78+0.64 


0.10 


27.9 


40.1 


35.2 


12.1 


32.7 


18.8 


24.489 


1.07E-05 




G58.78+0.65 


0.11 


27.9 


39.8 


35.2 


11.9 


32.5 


15.6 


20.091 


8.74E-06 




G59.79+0.63 


0.10 


27.5 


37.4 


35.5 


10.0 


34.5 


16.3 


24.380 


l.OOE-05 





^ References are for previously reported H2O masers with accurate positions from interferometric observations that fall within the polygonal 
aperture for the EGO published by''COl f ^3?3l . B02: Bcuthc r et~aI1 J2002D Bll: Bartkic wicz et al.| JgoTJ). BEll: IBreen & EUin gscn (mH^ CGll: 
ICaswell & GreenI | f20TB FC99: [Forstcr Caswell (.1999.) HC96: IHofner C hurchwellj i 19961) SlO: ISato et al.1 J2010D W06: [Wang ct al] J2006D 
indicates no H2O maser reference meeting these criteria was found in a SIMBAD search. 
"Outflow-only" source in |C08| . The published polygonal aperture does not include the "central" source, which is associated with H2O masers 
(G12.91-0.26; G34.26+0.15; G37. 55+0.20: IBreen & Eningsenll2011i : [Forster CasweIilfT999l : iBeuther eral|[200^ . respectively). 



Table 7 

H2O Maser Nondetections: Limits 



Source Name 




<J 


Liso 






Jy 






T 


Gil. 11-0. 11 


0, 


,14 


< 


1, 


.09E-06 


G12.02-0.21 


0, 


,10 


< 


1 


.41E-07 


G16.58-0.08 


0, 


,10 


< 


5 


.04E-08 


G16.61-0.24 


0, 


,10 


< 


5, 


.68E-08 


G19.36-0.03 


0. 


,16 


< 


3 


.40E-08 


G21. 24+0.19 


0. 


,10 


< 


1, 


.89E-08 


G24.11-0.17 


0. 


,09 


< 


9 


.87E-08 


G24.11-0.18 


0. 


,10 


< 


1, 


.09E-07 


G28.28-0.36 


0, 


,17 


< 


8, 


.07E-08 


G28.83-0.25 


0, 


,17 


< 


2, 


.OlE-07 


G28.85-0.23 


0, 


,12 


< 


1 


.73E-07 


G29.84-0.47 


0, 


,14 


< 


1 


.09E-07 


G29.89-0.77 


0, 


,16 


< 


1 


.74E-07 


G29.91-0.81 


0, 


,17 


< 


1 


.87E-07 


G29.96-0.79 


0, 


,06 


< 


6 


.52E-08 


G35.04-0.47 


0, 


,12 


< 


5 


.82E-08 


G35.83-0.20 


0, 


,13 


< 


2, 


.lOE-08 


G36. 01-0.20 


0. 


,10 


< 


1, 


.50E-07 


G40. 28-0.27 


0. 


,08 


< 


8, 


.94E-08 


G40.60-0.72 


0. 


,10 


< 


9, 


.26E-08 


G44.01-0.03 


0. 


,10 


< 


1, 


.16E-07 


G48.66-0.30 


0, 


16 


< 


4, 


,68E-08 


G49.27-0.32 


0, 


15 


< 


2, 


,18E-07 


G54.11-0.04 


0, 


16 


< 


1, 


,09E-07 


G54.11-0.05 


0, 


11 


< 


7, 


,35E-08 


G56.13+0.22 


0, 


,06 


< 


6, 


,80E-08 


G57.61+0.02 


0, 


11 


< 


1, 


,llE-07 


G58.09-0.34 


0, 


,06 


< 


1 


,69E-09 


G58.79+0.63 


0, 


,07 


< 


6 


,48E-08 


G62.70-0.51 


0, 


,06 


< 


4, 


.15E-08 



All limits are 4(t. 



Table 8 

Properties of Associated BGPS Dust Clumps^ 



Single Component^ Two Component'^ 



Source Name 


i:>CTjr S t.,al, iL) 


Si, 1,,, 


Distance* 


Radius 


Ra,(lius 


C'luiup ]Ma,ss 




C4uiii|) Mass 


T - 1 








(kpc) 


( ) 


I I'C) 




[Lo^^[<. Ill J J 


( Al . j 


l.l-oSl'i" I) 


GlO.29-0.13 


1474 


9.4(0.7) 


l..-i8(+0.86,-1.13) 


68.1 


0.52 


427 +^''* 

' —398 


4 1+'^'^ 

*--^-4.1 


585 


4.2 


GlO.34-0.14 


1483 


8.6(0.6) 


1.29(+0. 92,-1. 23) 


84.9 


0.43 


179 

^ ' ^ -179 


' -3.7 


361 


4.0 


Gil. 11-0. 11 


1589 


2.3(0.2) 


12.67(4-0.48,-0.41)'' 


56.7 


4.18 


11740 +^^*^ 

11 I'iV _2679 


3-ltlZ 


7967 


2.9 


G12.02-0.21 


1668 


1.0(0.1) 


5.30(-f0.20,-0.20)'5 


38.0 


1.75 


> 1024 till 


> 3.7tli 






G12.20-0.03 


1682 


2.1(0.2) 


11.70(H-0.31,-0.27)' 


52.8 


3.86 


-709 4-1674 


s.otll 






G12.42+0.50 


1710 


7.6(0.5) 


1.87(-|-0. 69,-0. 85) 


50.6 


0.62 


010 +384 
-239 


4 1+^-^ 






G12.68-0.18 


1747 


11.8(0.8) 


4.46(-|-0.30,-0.34) 


86.2 


1.47 


3454 tllf 


3.4li.? 


6447 


3.6 


G12.91-0.03 


1809 


3.6(0.3) 


4.48(-|-0.30,-0.34) 


65.0 


1.48 


1129 till 


3.213- 


1715 


3.4 


G12.91-0.26 


1810 


15.9(1.0) 


3.41(-|-0.43,-0.50) 


79.2 


1.13 


2615 tllll 


3-rtll 


3542 


3.8 


G14.63-0.58 


2081 


10.2(0.7) 


1.73(-|-0.63,-0.76) 


72.5 


0.57 


0(1 _42i 








G16.58-0.08 


2291 


2.6(0.3) 


3.24(-|-0.39,-0.44) 


77.6 


1.07 


927 tfrl 


3-3_3 1 






G16.59-0.05 


2292 


3.4(0.3) 


4.20(-l-0.30,-0.33) 


31.7 


1.39 


1148 ttll 


4.3lti 






G16.61-0.24 


2294 


0.8(0.1) 


3.39(-l-0.37,-0.42) 


28.7 


1.12 


195 tlf 


3.9it:? 






G17.96+0.08 


2358 


1.2(0.2) 


1.87(4-0.53,-0.61) 


26.4 


0.62 


. so +89 

> 82 


> 4.4ltl 






G18.67+0.03 


2431 


2.5(0.2) 


10.81(+0.24,-0.23)' 


56.7 


3.57 


4923 +'-2'35 


^■"'-2.3 






G18.89-0.47 


2467 


9.8(0.7) 


4.28(+0.28,-0.31) 


100.9 


1.41 


2266 till 


•'■"-2,6 






G19.01-0.03 


2499 


2.4(0.2) 


11.54(+0.29,-0.27))''''' 


44.1 


3.81 


4QQ7 +12S7 


3 1+^-^ 
•^■-^-2.6 


9056 


3.4 


G19.36-0.03 


2561 


6.8(0.5) 


2.09(+0.48,-0.54) 


95.7 


0.69 


438 till 


Q Q+3.8 

0.0_3 2 






G19.61-0.12 


2602 


3.1(0.2) 


3.84(4-0.31,-0.34) 


64.2 


1.27 


657 till 


9 + 3.1 
•^■^-2.9 






G19.61-0.14 


2602 


3.1(0.2) 


3.89(4-0.31,-0.34) 


64.2 


1.28 


645 till 


3.2l|i 






G19.88-0.53 


2636 


5.2(0.4) 


3.12(4-0.37,-0.41) 


40.7 


1.03 


718 tlt'r 


4.11^9 






G20.24+0.07 


2665 


1.2(0.1) 


4.38(4-0.28,-0.30) 


37.4 


1.45 


459 till 


3-etl.t 






G21.24+0.19 


2765 


1.1(0.1) 


1.92(4-0.47,-0.53) 


38.6 


0.63 


> 135 t\f 


> 4.lltS 






G22.04+0.22 


2837 


4.7(0.4) 


3.40(4-0.34,-0.37) 


87.9 


1.12 


728 tll'r 


3.0tl-° 


1376 


3.3 


G23. 01-0.41 


2971 


12.6(0.8) 


4.59(4-0.38,-0.33)' 


95.9 


1.52 


4024 tint 


, 9+3.1 

■J-^-2.9 


6401 


3.4 


G23.82+0.38 


3173 


0.6(0.1) 


4.43(+0.29,-0.30) 


18.6 


1.46 


> 84 ttl 


> 3.7ttl 






G23.96-0.il 


3202 


3.3(0.3) 


10.76(+0.28,-0.27)' 


57.3 


3.55 


5206 tllll 


2.9tl% 


14275 


3.4 


G24.00-0.10 


3208 


1.9(0.2) 


4.25(+0.29,-0.31) 


35.5 


1.40 


683 tfo° 


3-9tll 






G24.11-0.17 


3238 


2.5(0.2) 


4.62(+0.28,-0.29) 


77.4 


1.53 


1363 ttll 


3 0+^'^ 


1294 


3.0 


G24.11-0.18 


3238 


2.5(0.2) 


4.62(+0.28,-0.29) 


77.4 


1.53 


1330 ttll 


C).U_2 g 






G24. 17-0.02 


3246 


1.1(0.2) 


5.41(+0.29,-0.28) 


55.5 


1.79 


459 +1^0 

^■J" -146 


2.8tl:! 






G24.33-f0.14 


3284 


4.4(0.3) 


9.16(+0.29,-0.32)'' 


42.3 


3.03 




3.5tll 


4654 


3.5 


G24.63+0.15 


3383 


2.0(0.2) 


3.37(4-0.34,-0.36) 


50.0 


1.11 


641 1-4 


n 7 + 3.7 
•J- ' -3.4 


882 


3.8 


G24. 94-1-0. 07 


3440 


1.4(0.2) 


2.77(4-0.37,-0.40) 


49.3 


0.91 


126 til 


3.3_3 Q 






G25.27-0.43 


3481 


2.4(0.2) 


11.22(4-0.31,-0.30)*' 


66.8 


3.71 


> 9698 tllll 


> 2.9+11 






G25.38-0.15 


3503 


3.2(0.3) 


5.21(4-0.30,-0.29) 


33.5 


1.72 


1150 tfol 


3.9iri^ 






G27.97-0.47 


3862 


1.3(0.2) 


2.86(4-0.37,-0.39) 


60.1 


0.94 


301 till 


+ 3.5 

^■^-3.1 






G28.28-0.36 


3925 


6.3(0.4) 


3.05(4-0.36,-0.38) 


67.7 


1.01 


1478 till 


3-str, 


1458 




G28.83-0.25 


4055 


4.9(0.4) 


4.88(4-0.34,-0.33) 


64.9 


1.61 


1465 +^^^ 


3.2+3-1 






G28. 8.1-0. 23 


4056 


2.0(0.2) 


5.35(+0.38,-0.35) 


70.4 


1.77 


> 1615 +"'" 


> 3.otr, 






G29.84-0.47 


4230 


0.6(0.1) 


4.02(+0.34,-0.34) 


39.4 


1.33 


> 292 tl,f 


> 3.413- 






G29.91-0.81 


4245 


0.8(0.2) 


4.77(4-0.36,-0.34) 


27.1 


1.58 


> 731 +*3f2 


> 4.il|:^ 






G29.96-0.79 


4267 


2.5(0.3) 


4.81(+0.36,-0.35) 


61.5 


1.59 


1459 tltl 


3-3_3o 






G34.26+0.15 


5340 


78.5(4.8) 


3.52(+0.39,-0.39) 


103.4 


1.16 


9035 t^ 


3.813;^^ 


9123 


3.8 


G34. 28-1-0. 18 


5346 


0.8(0.1) 


3.38(4-0.38,-0.39) 


<16.5 


< 1.12 


189 tll° 


> 4.6ii:i 






G34.39-I-0.22 


5373 


20.8(1.3) 


1.56(4-0.12,-0.11)' 


96.1 


0.52 


001 +295 
°31 -225 


4.oirj 


1331 


4.2 



Table 8 — Continued 



Single Component 



Source Name 




o d 


Distanee^ 


Radius 


Radius 


C, 1 unip Mass 


Logfriff^] 






(3v) 




\ J 


(vc) 




G34.41+0.24 


5373 


20.8(1.3) 


1.56(+0. 12,-0. 11)' 


96.1 


0.52 


688 tf/e 


3-9tll 


G35.03+0.35 


5530 


6.5(0.4) 


3.24(+0.39,-0.39) 


89.0 


1.07 


814 till 


o 1 +3.2 
•^■-^-2.8 


G35.04-0.47 


5538 


4.7(0.4) 


3.11(+0.39,-0.40) 


98.5 


1.03 


> 1409 


. o + 3.4 
^ ■^••^-3.0 


G35.68-0.18 


5685 


3.4(0.3) 


1.76(+0.43,-0.45) 


87.5 


0.58 


350 _i8g 


o r + 4.1 
■J-0_3.4 


G35.79-0.17 


5700 


2.6(0.2) 


3.73(4-0.41,-0.40) 


54.3 


1.23 


476 +^"' 

-161 


3.3_3 Q 


G35.83-0.20 


5706 


1.3(0.1) 


1.80(4-0.43,-0.45) 


59.3 


0.59 


> 140 1^^" 




G36.01-0.20 


5722 


0.7(0.1) 


5. 48(4-0.58,-0. 58) 


36.1 


1.81 


> 436 till 


> 3.3tl:t 


G37.55+0.20 


5850 


3.6(0.3) 


5.67(-|-0.71,-0.71) 


51.2 


1.87 


qi)0O +1635 


q .7 + 3.8 

' -3.4 


G37.48-0.10 


5847 


0.8(0.1) 


3.59(-|-0.43,-0.42) 


37.0 


1.19 


> 222 tlT 


> 3.5_3 3 


G39.10+0.49 


5967 


0.8(0.1) 


1.47(-|-0.45,-0.47) 


49.5 


0.49 


> 29 +11 


> 3.5ttl 


G39.39-0.14 


5980 


1.4(0.1) 


4.22(4-0.56,-0.49) 


35.3 


1.39 


405 tlf. 


3.7+_ll 


G40.28-0.22 


6024 


3.5(0.3) 


4.98(4-0.64,-0.64) 


34.8 


1.64 


1008 


Q+4.0 

•^■^-3.6 


G40.28-0.27 


6023 


0.5(0.1) 


4.78(+0.59,-0.59) 


31.0 


1.58 


> 421 till 




G43. 04-0. 45(a) 


6110 


1.7(0.2) 


3.91(+0.66,-0.55) 


27.6 


1.29 


387 till 


4.iill 


G43.04-0.45(b) 


6110 


1.7(0.2) 


3.89(+0.65,-0.54) 


27.6 


1.28 


> 295 +\l\ 


> 3.9+_tl 


G44.01-0.03 


6134 


1.0(0.1) 


4.77(+0.77,-0.77) 


39.1 


1.58 


> 428 till 


> 3.4l|:^ 


G45.47+0.05 


6176 


5.2(0.4) 


4.77(+0.85,-0.85) 


46.5 


1.58 


1502 +'-"'^ 


3-rtti 


G45.47+0.13 


6177 


6.8(0.5) 


4.69(+0.82,-0.82) 


67.5 


1.55 


+1084 
J-'^J- -744 


3.3+11 


G45.50+0.12 


6177 


6.8(0.5) 


4.70(+0.83,-0.83) 


67.5 


1.55 


^ 5267 +3343 
^ 'f^y} 1 _2289 


> 3.8lti 


G45.80-0.36 


6202 


1.0(0.1) 


6.96(4-0.75,-0.75)'' 


25.7 


2.30 


> 632 ^229 


> 3.6+|-3 


G48.66-0.30 


6266 


0.9(0.2) 


2.41(4-0.59,-0.54) 


59.8 


0.80 


> 221 +111 


> 3A+11 


G49.07-0.33 


6298 


3.3(0.3) 


5.50(4-1.75,-1.75) 


63.9 


1.82 


1413 tl',f 


1+3.8 

^-^_3 


G49.27-0.34 


6323 


6.1(0.4) 


5.48(4-1.76,-1.76) 


51.8 


1.81 


2642 tllll 


s-eti.t 


G49.42+0.33 


6350 


0.8(0.1) 


12.02(4-0.56,-0.53)j 


27.5 


3.97 






G49.91+0.37 


6383 


0.6(0.1) 


0.53(4-0.52,-0.53) 


28.4 


0.18 


>etT 


> 4.8^ 


G53.92-0.07 


6444 


0.3(0.1) 


4.95(4-1.88,-1.88)j 


<16.5 


< 1.63 






G54.11-0.04 


6448 


5.3(0.5) 


3.73(4-1.00,-1.00) 


95.6 


1.23 


> 1750 tall' 


> 3.2+11 


G54.11-0.05 


6448 


5.3(0.5) 


3.69(4-0.98,-0.98) 


95.6 


1.22 


1316 tllf 


3.0tli 


G54.11-0.08 


6451 


7.0(0.6) 


3.85(4-1.06,-1.06) 


113.3 


1.27 


> 1287 t^fo" 


> 2.8ti:t 


G56. 13+0.22 


6464 


1.2(0.2) 


4.68(4-1.93,-1.93) 


44.8 


1.55 


> 988 tllf 


> 3.6ttl 


G57.61+0.02 


6472 


0.6(0.2) 


4.50(4-1.96,-1.96)' 


<16.5 


< 1.49 







Two Component'^ 
Clump Mass Log[n/f2] 
(M0) (Log[cm-3]) 



CO 
01 



Table 8 — Continued 



CO 
05 



Single Component^ Two Component*^ 
Source Name BGPS Cat. ID Si.imm^ Distance*^ Radius Radius Clump Mass Log[n/f2] Clump Mass Log[nj:f2] 
(Jy) (kpc) (") (pc) (Ma) (Log[cm-=^]) (Mg) (Log[cm-^]) 

^ Properties calculated as described in ^3.4| using nominal distances from Table[3] For EGOs matched to a BGPS source but undetected in NH3 emission, the BGPS clump ID and radius are 
listed, but no clump mass and density are calculated. Uncertainties in ^^'^ also in units of log(em~"'^). 

^ Calculated from the integrated flux density in the BGPS catalog assuming isothermal dust emission and T^ust —'^kin from Table [3] Quoted ranges include the uncertainty in the integrated 
flux from the BGPS catalog, the recommended BGPS flux correction factor, and the distance. See also ^3.4| and ^4.1.2| 

^ Calculated for sources fit with two NH3 components as described in ^3.4| Because the systematic uncertainties in this estimate are difficult to quantify, only nominal values are listed (which 
assume a BGPS correction factor of 1.5 and the nominal flux from the BGPS catalog). 

Flux density measurement taken directly from the BGPS catalog JRosolowsky et al.| j2010D . The recommended correction factor is applied when calculating clump masses and densities f ^3.4| | . 
^ This data is identical to that in Table [3] and is duplicated here for convenience. 

Associated with a 6.7 GHz CH3OH maser assigned the far distance bv I Green & McClure-Griffithi ) |201lD . Except for G12.20-0.03 and G45.80-0.36, all far distance assignments 
are "b" classifications in their scheme (see also i|3.2.1^ . We adopt the far kinematic distance estimated from the NH3 vlsr- 
^ The longitude and velocity of this source indicate that it is likely in the near 3 kpc arm (see for example Fig. 1 of I Green" 

iOnilSol)- Following [G reen McGlure-Griflathsl (i20TT|) . 
we place this source on a circle of radius 3.4 kpc around the Galactic Center, and adopt a distance uncertainty of ±0.2 kpc. 

It is unclear if this source is at the near or the far distance (e.g. ICvganowski et ani2011a| : IGreen fc McClure-Griffithsll20Tl) . 
' Masor parallax distance. References: G14.33-0.64. ISato et al.l ll20ia '). G23.01-0.41. IBrunthaler et af] ||2009(1 . G34.39+0.22 and G34.41+0.24: IKuravama et all ll2011f) . 
G35.20-0.74. IZhang et al.l 1120091) . 

NH3 nondetection. G49.42+0.33: distance estimated using H^^CO+ velocity from lC09l . G53.92-0.07: distance estimated using H2O maser peak velocity. G57.61+0.02: distance ^> 
estimated from velocity of weak {3.9(t) NH3(1,1) emission below our 4cr detection threshold. The distance for G57.61+0.02 is included here and in Figure[^for completeness, but ^ 

this source is otherwise excluded from our analysis. See also ti3.2.1l '2, 

k O 

The density range is not constrained because the lower end of the distance range is 0.0 kpc. ^ 

s 

H 
> 



Table 9 

H2O Mascr and NH3 Surveys of MYSO Samples: Detection Rate Comparison 



Sample 


Reference^ 


Sensitivity 


H2O Maser 
Resolution*^ 


Detection Rate 


Sensitivity 


NH3(1,1) 
Resolution^ 


Detection Rate 


EGOs (overall)"^ 
EGOs (IRDC)'' 
r!ji^ ws i^v^iass i K^n.-^\jri) 
EGOs (Class II CHsOH)'^ 
EGOs (Class I & II CH30H)*= 
EGOs (Neither CHsOH)"^ 


this work 


a ~ 110 mjy 


SD (73") 


68% 
62% 
yu /o 
86% 
95% 
33% 


(7 ~ 50 mK 


SD (73") 


97% 
100% 
1 nn'^^^ 

iUU /o 

96% 
100% 
93% 


Dust Clumps (overall)*^ 

Uust Clumps (mm only) 

Dust Clumps (Class II CH3OH))" 

Dust Clumps (Class II CH3OH & cm cont.))'' 

Dust Clumps (cm cont.))" 


BEll/HlO 
BEll/HlO 
BEll/HlO 
BEll/HlO 
BEll/HlO 


a ~30-40 mJy 


Int. (~2") 


44% 
23% 
79% 

83% 

58% 


variable 


SD (52") 


56% 
53% 
65% 
50% 
67% 


Dust Clumps (overall) 


Dll 


(T ~ 60 mJy 


SD (33") 


40% 


a ~ 100 mK 


SD (30") 


72% 


RMS (North) 
RMS YSOs (North) 
RMS HII (North) 
RMS YSO/HII (North) 


UU 
Ull 
Ull 
Ull 


cr ~ 120 mJy 


SD (33") 


52% 
52% 
52% 
56% 


(T ~ 50 mK 


SD (30") 


81% 

85% 
78% 
85% 


RMS (South) 
RMS YSOs (South) 
RMS HII (South) 
RMS YSO/HII (South) 


U09 
U09 
U09 
U09 


cr ~ 250 mJy 


SD (2') 


27% 
26% 
28% 
32% 








IRDC cores (overall) 
IRDC cores (quiescent) 
IRDC cores (active) 
IRDC cores (red) 


Ch09 
Ch09 
Ch09 
Ch09 


(T ~ 40 mJy 


SD (33") 


35% 
16% 
59% 
54% 








ffl^S-select HMPOs 


Sr02 


a ~400 mJy 


SD (40") 


42% 




SD (40") 


86% 


fflylS-selected UC HIIs 


C90 


detect, limit 400 mJy 


SD (40") 


67% 


detect, limit. 300 mK 


SD (40") 


70% 


H2O masers 


A96 








variable 


SD (1.4') 


58% 


OH masers 


BCEPIO 


CT ~40-150 mJy 


Int. (~2") 


79% 








Class II CH3OH masers (6.7 GHz) 
Class II CH3OH masers (6.7 GHz) 
Class II CH3OH masers (6.7 GHz) 
Class II CH3OH masers (6.7 GHz) 


Bll 
Sz05 
B02 
P12 


cr ~2-10 mJy 
a ~450 mJy 
a ~ 1000 mJy 


Int. (~0'.'15) 
SD (40") 
Int. (~1") 


71% 
52% 
62% 


CT ~ 30 mK 


SD (40") 


95% 



^ A96: lAnglada et al.Nl996r . B02: iBeuther et al.l (f200^ . BCEPIO: fBreen et all (l2010bD . Bll: iBartkiewicz et anil201lD . BEll: [Breen & EllingsenI ll201ir . Ch09: [Chambers et a'nil200gD . C90: 
IChurchwell et all (fT990h . Dll: I Dunham eTall i2011br) . HIO: iHill et all ||2010D P12: iPandian eraP <f201^ Sr02: ISridharan et al.l ll2002f) . Sz05: ISzvmczak et all i2005ll . U09: lUrauhart et all 
| f200g|) . Ull: lUrquhart et all J2011D 

^ Int. indicates interferometric observations, and estimated positional accuracy is given. SD indicates single dish observations, and the FWHP beam size is given. 

^ See Table [2] for additional statistics on EGO subsamples. 'Class I': all EGOs with a Class I maser detection in ICElH , regardless of Class II association (or lack of Class II information). 
'Class II': all EGOs listed as Class II maser detections in lCElll . regardless of Class I association. 'Class I & II CH3OH': EGOs listed as both Class I and Class II detections in lCElH . 'Neither 
CH3OH': EGOs listed as both Class I nondetections and Class II nondctections in lCElH . 

BEll: H2O maser survey; HIO: NH3 survey. These surveys both target subsets of the iHill et al.l l |2005r dust clump sample, and divide their observed dust clumps into the same general 
categories. The studies are, however, independent, and BEll do not include the NH3 results of HIO in their analysis. 



CO 



